NOT TO BE TAKEN FROM THIS ROOM 


oe Ota ee AO" = aoe’: 


ae ret oee tes oe os 








Gx apnis 
UNITASITAIS 
LARERTAEASIS 





COO 


Digitized by the Internet Archive 
In 2021 with funding trom 
University of Alberta Libraries 


https://archive.org/details/Fuhr1974 








oO 9 Aad 8S & @ 





WAME OF Alrriioe, anya Jou n 
TITLE OF THESIS: WUCLZAR MACEPCIC aa 


OF THe Six W 
DEGREE FOR WHICH THESTR war, 
FEAR THIS COGHEE Saeed: 1474 


Perniasion is taraby 
SENTA LIHRARY ¢ repend 
thesia wand to ten’! w ‘) ow 


ei. 


schooler iy is Ssciantbi 


The author remem: otmer par Winn cight a 
heither the thesia one - ei pom it 
may Se printed or oth . oF 


author's #wri*®ten pecei-e 








1 kas WANS yey Ven othe ey Tony Oi Bot Bon Rf A 


RELEASE FORM 


NAME OF AUTHOR: BRYAN JOHN FUHR 

TITIES OFOTEESIS : NUCLEAR MAGNETIC RESONANCE STUDIES 
OF THE SOLUTION CHEMISTRY OF METAL 
COMPLEXES 

DEGREE FOR WHICH THESIS WAS PRESENTED: DOCTOR OF 
PHILOSOPHY 


YEAR THIS DEGREE GRANTED: 1974 


Permission is hereby granted to THE UNIVERSITY OF 
ALBERTA LIBRARY to reproduce single copies of this 
thesis and to lend or sell such copies for private, 
scholarly or scientific research purposes only. 

The author reserves other publication rights, and 
neither the thesis nor extensive abstracts from it 
may be printed or otherwise reproduced without the 


author's written permission. 


i 


ataueauva OO FRETS 





aug veew WAYRE ‘SOMO TO AMAM 

oarquv2 SAMOS DTTIMOAM AARIIGM :BYSHn? IO WeTTT 
jatam To yarTerMand KOTTUUCe aut TO 
2K aad IM0D 

2o ROTDOd ;aeTaesits saw SfaaRe HOLHW #OT audoed 


YROOSOITNS 


6fol Gata saneed arnt #AaY 


qo vera ayryu UAT. ot, bagditarp yueu2s8 2i aolee Lined 
Hint 2O e9iyoo Sibare, souboImges o2 YAATSIA AHIGIR 
stnyita tok 2atgoo doue {i16@ xo Bael of Das ateed? 
vine esanqiie dptesres otitdieive wo ylusiodce 


eintpis woissoridgg 7 five p38S2 sonsve oft 








THE UNIVERSITY OF ALBERTA 


NUCLEAR MAGNETIC RESONANCE STUDIES OF 


THE SOLUTION CHEMISTRY OF METAL COMPLEXES 


BY 


(C) BRYAN JOHN FUHR 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES 
AND RESEARCH IN PARTIAL FULFILMENT OF 
THE REQUIREMENTS FOR THE DEGREE OF 


DOCTOR OF PHILOSOPHY 


DEPARTMENT OF CHEMISTRY 


EDMONTON, ALBERTA 


SPRING, 1974 


Sv 
1 7 e 


ATasata FO VEVAR 


19 aa TOUTS OVA ORM arenas 





(4fEMOD LATAM Wh YATSa 


% ' : | 


1T2eaHyT A 


f | , 4 cy? 
o - + & 4 ~/ - 
e ge "y ar re < aT 
ri, - u q 4 Lida . {J on 5 
» m en 3 2% ae 9 
225 4 « 4 fa Ad Find te s be 4 a 
' ‘ ; 
& ) tHe * \ J A) — - 
—=~ 
7 a 
¥ mi *seeoU0r 7 ane . every a] 
ivi iGk 2 4 i) i” IMTHAN: () i 
4 ran a 





THE UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, 
and recommend to the Faculty of Graduate Studies 
and Research, for acceptance, a thesis entitled 
NUCLEAR MAGNETIC RESONANCE STUDIES OF THE SOLUTION 
CHEMISTRY OF METAL COMPLEXES 
submitted by BRYAN JOHN FUHR 
in partial fulfilment of the requirements for the 


degree of DOCTOR OF PHILOSOPHY. 


es 43) 0 a 






avanara 0 ¥ | 


beer sved yes sada yisoeee poupiasghay sat 
asibyse gtaubextd 20 yi ives? eds o¢ Sasamooes brs 
, sored eeest is} Aores2e8 






bhefsiens abgens 5 













7 |. tae 2 lee 


ams « Part t 


aq": - 
Proton magnetic resonance iper) methods have 
. 4 _ 
been used tr teri ne the evartnge 4 7 et haw ’ J the 
td 
nonequivalent SGetate® Gethyienic protons betore 211 - 
_ 


enange in the cadmiun ané eine comp! ce: <J- 
propylienediaminetetraacetic aci« , 2-PDTA)  &n 
t,é-baty! enedianinetetraacetica acid (1,4~B80TA 
Srgecisvelchce a3 indicated by AS mul! 

im the pmr opectra of these com Kis 
Were Geternined ever a nae 7 

the extent of opliarse of. the AU puerttern: yy 
experimenta! ‘Speetrs with theore mal ebectra 


were calcujated ag a fancteson of the avers cde 
* 


TO MY WIFE, ARLENE 
before An inter thange. a 2 ie a 


e Se 


bifetines are 4 measure of the rates of partie! 


Gissociation of these complexes, The AB multip 
i 


patterns in the per aepectra of the cudmium, sinc 
afd lead ctuplexos of ethy lenediesine-N-N' <diacetic- 
Nel'-Gigropionic acid (EDDDA) were also invesatigqeted 
aoa) function of tempersture. ms 


a . - : 7 
: fh of the liga 7 4Cka : 
. a ee me : _ eS 7 ' ; eae oe de} - 


> 
9 ra annids 
’ . a 


a an! fe ae) 


oe ean i) a 
: Rae 


“—— 


- 





ABSTRACT 


Part I 


Proton magnetic resonance (pmr) methods have 
been used to determine the average lifetimes of the 
nonequivalent acetate methylenic protons before inter- 
change in the cadmium and zinc complexes of 1,3- 
propylenediaminetetraacetic acid (1,3-PDTA) and of 
1,4-butylenediaminetetraacetic acid (1,4-BDTA). 
Nonequivalence is indicated by AB multiplet patterns 
in the pmr spectra of these complexes. The lifetimes 
were determined over a range of temperatures from 
the extent of collapse of the AB patterns by matching 
experimental spectra with theoretical spectra which 
were calculated as a function of the average lifetime 
before AB interchange. It is proposed that these 
lifetimes are a measure of the rates of partial 
dissociation of these complexes. The AB multiplet 
patterns in the pmr spectra of the cadmium, Zinc 
and lead complexes of ethylenediamine-N-N'-diacetic- 
N-N'-dipropionic acid (EDDDA) were also investigated 
as a function of temperature. 

The kinetics of the ligand exchange reactions 
of the cadmium, zinc, and lead complexes of 1,3-PDTA 


and EDDDA have been studied using pmr line-broadening 
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methods. Rate constants were obtained for the proton- 
assisted dissociation of each of these complexes 

and for the displacement of complexed ligand by free 
ligand for several of the complexes. The rates of 
first order dissociation of all the complexes are 

too slow to measure by the pmr technique. Rate 
constants for the reaction of cadmium and zinc 

with monoprotonated 1,3-PDTA are approximately 

equal to those of the analogous reactions with 
monoprotonated ethylenediaminetetraacetic acid, 

while the rate constants for reaction of cadmiun, 
Zinc and lead with monoprotonated EDDDA are from 

2.5 to 10 times slower than the rate constants for 
reaction with monoprotonated 1,3-PDTA. Mechanisms 
are proposed for the formation reactions of the mono- 


protonated ligands. 


Parco. 


The binding of zinc, cadmium, lead and mercury by 
the tripeptide glutathione has been investigated 
by oc magnetic resonance spectroscopy. Binding to 
the potential coordination sites was monitored as 
a function of solution conditions by observing the 
chemical shifts of the carbon atoms of glutathione. 


The results indicate that each of these metal ions 
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binds to the potential coordination sites of gluta- 
thione with a high degree of specificity. Mercury 
binds only to the sulfhydryl group at mercury to 
glutathione ratios up to 0.5 and over the entire 
accessible pD range. At a metal to glutathione 
Pastor) 0.5;) 21nC and caumium bind to both the 
sulfhydryl and amino groups, the extent of binding 
to the two different sites being a function of pD, 
while lead binds only to the sulfhydryl group. 

Some binding of the carboxylic acid groups to 

Zinc, cadmium and lead was detected in certain pD 
regions. The chemical shift data also suggest that 
Zinc-promoted ionization of the peptide protons 
with subsequent binding of zinc to the ionized 
peptide nitrogen takes place above pD 10.5. The 
results are discussed in terms of the possible 


structures of the complexes. 
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CHAPTER TI 


INTRODUCTION 


A. DISSOCIATION KINETICS 


Formation and dissociation reactions involving 
multidentate ligands proceed via stepwise mechanisms 
in which one of the steps governs the overall rate 
of the reaction (1). For example, Ahmed and Wilkins 
(2) found that che Pioneeic results for the dissocia- 
tion of the nickel complex of ethylenediamine were 
consistent with a mechanism involving stepwise rupture 
of the individual dentate bonds. Margerum and co- 
workers (3), on the basis of kinetic data for the 
formation of some nickel-polyamine complexes, proposed 
a mechanism involving the stepwise formation of the 
individual coordinate bonds and concluded that the 
formation of the first nickel-nitrogen bond was the 
rate-determining step. It was possible for these 
workers to estimate the values for the stepwise 
equilibrium constants from their kinetic data. 

Stepwise mechanisms have also been proposed for 
reactions involving the exchange of a multidentate 
ligand between two different metal ions and of a 
metal ion between two multidentate ligands (4). For 


instance, such mechanisms have been inferred from 
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kinetic data for the reactions of zinc and copper 
with the nickel complex of ethylenediaminetetraacetic 
acid (EDTA) and analogs of EDTA (5,6,12,19) and for 
the reactions of EDTA with several nickel-polyamine 
complexes (7,8). In particular, the proposed mechanism 
for the reaction of copper with the nickel complex 

of EDTA (6) involves the unwrapping: of an iminodi- 
acetate fragment from the complex, followed by reaction 
with copper to give a dinuclear intermediate. For 
this case, the kinetic data yielded a rate constant 
for the half-unwrapping of EDTA from nickel. 

The preceding work, in which such techniques 
as isotopic labelling, spectrophotometry and tempera- 
ture jump relaxation were used, has shown that exchange 
reactions involving metal complexes of multidentate 
ligands take place in a stepwise fashion. However, 
it has only been possible in a few cases, such as 
in the work of Margerum and coworkers (6), to measure 
the rate constants for the individual steps. 

Proton magnetic resonance (pmr) studies of metal 
complexes of EDTA) (9-11,13-17) and analogs of EDTA 
(10,13,15-18, 20-23) have shown that the complexity 
of the pmr spectrum of complexes of this type 
depends on the lability of the individual metal- 


ligand bonds. Hence pmr may provide a way to study 
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the kinetics of individual metal-ligand bonds. EDTA 
contains both oxygen and nitrogen donor groups and, 
depending on the lability of the various metal-ligand 
bonds, several situations may arise. The first is 
the one in which both the metal-oxygen and metal- 
nitrogen bonds are labile on the pmr time scale. In 
this case, because of rapid nitrogen inversion, the 
acetate methylenic protons are all equivalent, and 
the pmr spectrum for these protons consists of a single 
resonance peak. An example of a complex giving rise 
to this type of spectrum is Sr (EDTA)“ (39). 

The second case, in which the metal-oxygen bonds 
are labile but the metal-nitrogen bonds are inert on 


the pmr time scale, may be illustrated by the 


following rotational conformers. 


Et Ney ES Ac Et Ac 


on praca <e 


: COmmicos H 
M 


The acetate group is free to rotate about the nitrogen- 
acetate methylenic carbon bond. Examination reveals that, 
regardless of the rate of rotation, leading to inter- 


change of the above three conformations, Hy and Hy do 
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not necessarily experience equivalent environments 
because of the asymmetry of the quaternary nitrogen 
atom. Hence Hy and A, may have different chemical 
shifts and if so, their pmr spectrum will consist of 
a single AB pattern. Interchange of Hy and H, can 
occur only by inversion of the nitrogen atom, which 
presumably is possible only when the nitrogen atom is 
not metal coordinated. As long as the rate of inter- 
change of Hy and Hp is less than the chemical shift 


difference Adan (in sec 1) between H. and H an AB 


A B’ 
pattern will be observed. In other words, metal- 
nitrogen bonding is inert on the pmr time scale. The 
pmr spectrum of the acetate methylenic protons of 
Cd (EDTA) 77 consists of a single AB pattern such as 
that described above (24). 

For the situation in which both the metal-oxygen 
and metal-nitrogen bonding is inert on the pmr time 
scale, rather complex spectra are expected because 
the acetate groups are no longer equivalent. In the 
case of an octahedral complex in which all the ligand 
atoms are bound, two of the acetate groups are in 
the ee formed by the metal and the two nitrogen 
atoms and two are out of the plane. Therefore the 
acetate groups may be expected to exhibit two different 


AB patterns in their pmr spectra. The inert bonding 
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in Co (EDTA) 77 results in two different AB patterns 
in the pmr spectrum for the acetate methylenic protons 
Gf this -coanp bax "(9) ; 

Sudmeier and Reilley (24) obtained evidence for 
a stepwise mechanism for the reaction of EDTA with 
Ca (EDTA) 77 by monitoring the degree of collapse of 
the AB pattern in the pmr spectrum of the complex 
under changing solution conditions. In particular 
they studied the relative effectiveness of ligands 
which simulate fragments of EDTA in causing AB 
interchange. At some point in the overall displace- 
ment reaction, the first nitrogen of the complexed 
EDTA ligand is displaced. Dissociation of this metal- 
nitrogen bond was studied from the collapse of the 
AB pattern. It was found that ligands most effective 
in causing AB interchange contain at least three 
coordinating groups, a nitrogen and two carboxylate 
groups. On this basis, it was postulated that formation 
of the third bond is most effective in breaking the 
metal-nitrogen bond of the original complex. A 
stepwise mechanism was proposed with the rate-deter- 
mining dice being the breaking of the first metal- 
nitrogen bond of the originally complexed ligand and 
displacement by a carboxylate group of the incoming 


ligand. 
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Information about the kinetics of dissociation 
of the proton from the nitrogen in acidic solutions 


of tertiary amines of the structure 


hie 
(ip ata 
ane Mina 

He, = 


has been obtained from collapse of the AB pattern for 
the methylenic protons (25-28). When the nitrogen 
in these types of amines is protonated, it cannot 
undergo inversion. In this case, the following three 
rotational conformers are important, assuming rotation 


about the carbon-nitrogen bond. 


As in the case of the EDTA type complexes, when the 
rate of nitrogen inversion is slow on the pmr time 
scale, H, and Hp will not necessarily experience 
equivalent environments regardless of the rate of 


interconversion of the above three conformers. Hence 


the pmr spectrum for the methylenic protons will 
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consist of an AB pattern. As the inversion rate is 
increased by increasing the pH of the solution, thus 
reducing the fraction of protonated amine in solution, 
the rate of interchange of Hy and H, is increased 
resulting in a collapse of the AB pattern to a single 
peak. The rate of this AB interchange was determined 
for a series of tertiary amines by lineshape analysis 
of the pmr spectra (27). It was shown that if the 
rate of nitrogen inversion is very much less than the 
rate of reprotonation of a deprotonated amine, the 
observed AB interchange rate will be some fraction of 
the rate of nitrogen inversion. However, if the rate 
of inversion is very much greater than the rate of 
reprotonation, the observed AB interchange rate will 
be one-half the rate of proton exchange. 

The work of Sudmeier and Reilley (24) and that 
on the tertiary amines (27,28) indicates that from a 
lineshape analysis of AB patterns it may be possible 
to obtain detailed information about the dissociation 
kinetics of complexes of multidentate ligands. By 
analogy with the work on tertiary amines, collapse 
of the AB’ pattern by inversion of a nitrogen atom in 
a complexed aminocarboxylic acid ligand may be governed 
by the rate of metal-nitrogen bond dissociation and/or 


the rate of inversion of the noncoordinated nitrogen 
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atom. If dissociation of the metal-nitrogen bond 
governs the rate of collapse, the presence of an AB 
pattern is an indication of inert metal-nitrogen 
bonding and potentially can provide information about 
the partial dissociation kinetics of the complex. 

In Chapter II, the results of a pmr temperature 
study of the AB multiplet patterns arising from the 
acetate methylenic protons of the cadmium and Zinc 
complexes of 1,3-propylenediaminetetraacetic acid 


(1,3-PDTA), 


HO. CCH CH.COlE 
NCH Co. CHAN 


Ya 
iia 
eee 


HO,CCH., CH,CO,H 


Mya ated Ws 


and of 1,4-butylenediaminetetraacetic acid (1,4-BDTA) 


are presented and discussed (29,30). These particular 
aca NS CH.CO,H 
ra NCH,CH,CH,CH.N 
HO,CCH, CH,CO,H 
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ligands were chosen since they are structurally 

Similar to EDTA whose cadmium complex, as noted above 
(24), exhibits an AB pattern in its pmr spectrum. 

It was thought that the different sizes of the nitrogen- 
metal-nitrogen chelate rings in the metal complexes 
OfCL,3-PDTA and’ 1,4=-BDTA might affect the lability 

of the metal-nitrogen bonds. The mutliplet patterns 

for the acetate methylenic protons of the cadmium, 

zinc and lead complexes of ethylenediamine-N,N'- 
diacetic-N,N'-dipropionic acid (EDDDA) were also studied 
as a function of temperature. 


HOCCH CH, CH,CH,CO,H 


NCH.CH,N 


HO,CCH. CH,CO.H 


EDDDA 


The main objective of this work was to determine what 
kinetic information can be obtained from the multi- 
plet patterns of the methylenic protons in complexes 


of this type. 


B. FORMATION KINETICS 


According to the dissociative mechanism for the 
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reaction of an aquated metal ion with a ligand to 

form a metal complex (31-33,70), the metal ion and the 
ligand first diffuse together to form an outer-sphere 
complex. A water molecule then dissociates from the 
aquated metal ion and a donor groupv of the ligand 
coordinates with the metal ion at the vacated site. 
When steric and electrostatic effects are negligible, 
the rate-determining step in this reaction mechanism 
is water loss from the aquated metal ion. For a 


monodentate ligand this mechanism may be depicted as 


k 


M(H ore = 3h USPeMiH 0) 27 . Smee 
aE Ks) ae tes 
-x 
K 7 H20 - 
> M(H,0) .X + H,O 


For a multidentate ligand, this reaction scheme has 

to be modified to account for the stepwise nature 

of complex formation with multidentate ligands (2,7,8). 
This is illustrated by the reaction scheme for a 


bidentate ligand. 
H50 
wae. ie soe 
le > 2S ge Coyie Sivee ts exe 
M(H,O)¢ + | 2 M(H,0)¢ X—X <> H50 
x Ky 


k 
2 
M(H.O) -X—X os M(H.,O) , (X—X) 
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As in the case of a monodentate ligand, the bidentate 
ligand forms an outer-sphere complex with the metal 
ion, followed by water loss from the metal ion and 
formation of the first metal-ligand bond. A second 
water molecule then dissociates from the metal ion in 
a position cis to the first metal-ligand bond, and 
the second metal-ligand bond forms. The rate of this 
second bond formation process is represented by ky in 


the above reaction sequence. Depending upon the 


H50 
i 


the first or second bond to the bidentate ligand 


relative magnitudes -of k : and Koy the formation of 
determines the overall rate of complex formation. For 
multidentate ligands such as EDTA, the reaction se- 
quence for complex formation would include steps for 
each successive bond formation. 

It has been found for the majority of complex 
formation reactions involving multidentate ligands, 
that the formation of the first metal-ligand bond is 
rate-determining (70), its rate being governed by 
the rate of water dissociation from the aquated metal 
ion in the outer-sphere complex. For example, the 
rate of eon eee on DE. divalent cadmium, 2inc, 
cobalt, copper and nickel by the monoprotonated form 
GietDUA CLS, 24,54, 36-55). | 


- - + 
m?* 4 wEpTA? > M (EDTA) * + H 
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are all approximately equal to those predicted by 
the dissociative mechanism. 

Other studies have shown, however, that the rates 
of complexation of calcium and strontium by the 


tetraanion of EDTA are factors of 2.5 x 102 and 


Me ahaBDTR ac M (EDTA) 77 


See ibe greater than their rates of complexation by 
HEDTA>~ (39740). Similarity, the rates of complexation 
of nickel (41) and cadmium (42,44) by the trianion of 


NnLtrilotriacetic acid (NTA) are 7 x 104 times 


Mou LEANDRA tas M(NTA) _ 


greater than their rates with HNTA*”, while the rates 


Of reaction jot zinc (44) and lead (43) with HNTA?~ 
are three to four orders of magnitude less than 

the approximate values predicted by the dissociative 
mechanism. 

In one of the mechanisms proposed to account for 
these decreased rates, a small amount of the mono- 
protonated ligand is considered to be present in the 
carboxylate-protonated form along with the nitrogen- 
protonated form. This carboxylate-protonated form 
is considered to be the reactive species due to the 


Blocking of the nitrogen atom by the proton (44-46). 
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A second reaction mechanism to account for the decreased 
rates of reaction of Meas with om for ligands of the 
EDTA type may be illustrated by reference to Figure l. 
Reaction proceeds by dissociation of water molecules 
followed by metal-ligand bond formation to yield 
intermediate (III) in which the three dentates of one 
end of the ligand are metal coordinated while co- 
ordination of the other end is blocked by the proton 
OneLhe nitrogen,  —roryiurther coordination totoccur, 
the proton must be transferred to a carboxylate 
oxygen or a solvent molecule before the intermediate 
dissociates to the reactants. A slow proton transfer 
step relative to dissociation of the intermediate 
back to reactants would account for the decreased 
rates mentioned above. 

The size of a chelate ring formed in many 
complexation reactions may also affect the overall 
rate of the reaction..ivor instance, the rate of 
reaction of cobalt(II) with o-alanine, in which a 
five-membered ring is formed, is an order of magnitude 
greater than that for the similar reaction with 
ee aine in which a six-membered ring is formed (47). 
Similarly, the rates for the reactions of nickel and 


cobalt with a-aminobutyric acid were found to be 


about an order of magnitude greater than the corresponding 
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rates with the g-amino acid (48). The results with 

the a-amino acids are consistent with the dissociative 
mechanism. With the f-amino acids, however, the results 
are consistent with a mechanism in which ring closure 

is the rate-determining step. The latter mechanism 
implies that the repossession of the vacated coordination 
site by a water molecule is a faster process than ring 
closure. 

The size of chelate rings formed may also be an 
important factor governing the rates of the formation 
reactions of EDTA type complexes from the mono- 
protonated form of the ligand. For example, the larger 
the chelate ring which must be formed in going from 
intermediates (IV) to (V) in Figure 1, the more 
likely the vacated coordination site is to be reoccupied 
by a water molecule before ring closure can occur. 

The stability of certain intermediates may also affect 
the overall rates of these formation reactions. For 
instance, the less stable intermediate (III) is, the 
more likely 1t 18 to dissociate to reactants before 
proton transfer and bonding of the second nitrogen 

can occur. In order to see if these two factors do 
have an effect on the rates of certain complexation 
reactions, a pmr line-broadening study of the kinetics 


of the ligand exchange reactions of the cadmium, zinc 
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and lead complexes of 1,3-PDTA and EDDDA was under- 
taken (49). The results of this study are described 

in Chapter III. The predominant pathways by which 

the ligands exchange between free and complexed forms 
are established and their rate constants measured. 

The results are discussed with emphasis on the mechanism 
of formation of the complexes from the monoprotonated 


forms of the ligands. 
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CHAPTER II 
PARTIAL DISSOCIATION KINETICS OF METAL COMPLEXES 


OF AMINOCARBOXYLIC ACIDS 


in this ichapter; the results of a study“of “the 
temperature dependence of the AB multiplet patterns 
in the pmr spectra of certain metal complexes of 
aminocarboxylic acids are presented and discussed. 
In particular, the cadmium and zinc complexes of 
both 1,3-PDTA and 1,4-BDTA and the cadmium, zinc and 
lead complexes of EDDDA were investigated. A central 
purpose of this work was to ascertain what kinetic 
information can be obtained from the nonequivalence 


of the acetate methylenic protons in complexes of this 


type. 


As eC RESULTS 


Analysis. of Multiplet Patterns 

The 60 MHz pmr spectrum of the acetate methylenic 
protons of zn (1,3-PDTA) 7 in aqueous solution at 25° 
is shown, in the lower left hand side of Figure 2. 
The pH of the solugion was %8.7 at 25).% The spectrum 
consists of an AB pattern, the upfield peak of which 
is obscured by a broad triplet due to the two 


methylene groups bonded to the nitrogen atoms in the 
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' propylene part of the ligand. Analysis of this AB 
pattern by standard methods (35) yielded the values 
for Asan and Ing listed in Table I. Asan is the 
chemical shift difference between the A and B protons 
and Jap is the spin-spin coupling constant. The 
value for this geminal coupling constant was assumed 
to be negative (17). 

The 100 MHz pmr spectrum of the acetate methylenic 
protons of Ca (1,3-PDTA) 7 in D,0 solution at 10° is 
shown in Figure 3. The pD of the solution was 7.5 
at 25°. This multiplet pattern is more complicated 


than that for the zinc system because of the eae 


and ree isotopes which have a nuclear spin quantum 


d 


number of 1/2 and are present at a combined natural 
abundance of 25.0%. The multiplet consists of the 
AB part of an ABX pattern superimposed on an AB 
pattern. The AB pattern arises from the complexes 

of the cadmium isotopes whose spin numbers are zero, 
while the AB part of an ABX pattern is due to the 
cd and ed complexes. The broad triplet just 
upfield from the multiplet pattern arises from the 
two methylene groups bonded to the nitrogen atoms 

in the propylene part of the ligand. Analysis of the 


AB pattern from slow exchange spectra of the type shown 


in Figure 3 at several temperatures in the range 





























7 
- 


oe ¥ 
aA eids to abeylenk vbampe nat oF 
ah an hel 
ouiny, 984 bebteiy (¢£) aboitsau i stinied yd nett 
_ 
sity at _,bA «1 @hdat mr. podekt.. a 
Sw aa” ® i 
non 
if caowlad soeata pride fse net C 


zaazgo sg 8 Gre A 96) 
4?  .4nazenoo onitauoo Snes: Pers at 
b saeran caifawos feminap aids Pr’ ‘wom 
AVE) ort rnEes od 
ps 
jo ott Yo mvatoogs Say SAM OGL = ait 7 
of 2 ui az * (ANAGHE . TYBD Jo: enosone 
nettirtoe et7 tote ait 6. Saye w sie 
) nettmbtoe si : I t ie i stoi 
Y o = 7. t os * F i vn p 7 
Son el aryetv6q Joteqss.an = cit ‘ es. - 
, - 
nde 26 genncesd mateve belts ons ~ot tsadd fee 
Z _ 
: iat ee 
mjnpve arga «xsefouvn 6 Svs NOE a9co7oOeL bo ~ fine 
; i, i ‘oa 
nidmods ns +6 anges sis bis <\0 te seca 
! shia 
”~ + Ma _ oe : 7 
Sils shop seigisium sHt . .f0.eS 20 sonebowds 


pocoumiyeque musiteq ASA 20 srsq G6 
s+ mon? eoebas adedzag GA ect xesis 
SaqosogL minhes ont to 


'} ot onb ab cistiaq, X6A as 20. 32eq A orl, =! 
ELL tf 


saut ssfuqizt BPeotd off -epxelqmon BS Bris 
a 7 
f+ wert es@itea niset¢eg sotaldium ais mos? Biakaae 


amod 6 naporsin satis oj bebaod a one etoly 

_ Ve Lhe 

a Wiel o bentacaie _ is 26 i ams wyTO 

= (ee Bae shaw M : ata Ly qosg 
ae i” ae ( 

















Table I 
Pmr Spectral Parameters for the Acetate Methylenic 


Protons of the Indicated Metal Complexes 


ds,p(ppm) Jy, (Hz) IF ax! 
Cd (1,3-PDTA) 77 0.192 -16.3 1128 
gn(1,3-PDTA) ~~ 0.220 -16.3 
Pb (1, 3-PDTA) > i 19.2 
Ca (1,4-BpTa) 77 Or141 -16.2 1728 
ee ene 0.204 S152 
Cd (EDDDA) 77 b b 13.5 
zn (EDDDA) ~~ 0.453 -18.3 
Pb (EDDDA) 77 Oee21 -16.7 Part 


(a) From reference 30. 


(b) Could not be determined from the spectrum. 
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3 to 26° yielded the value for J listed in Table I. 


AB 
The chemical shift difference was found to decrease 
linearly as the temperature was increased, and a 


least squares analysis of the data yielded the 


following relation: 


45. , (ppm) = -5.00 x ee Pee. 200 I 


where T is the temperature in degrees centigrade. 
Analysis of the AB part of the ABX pattern using the 
value determined for Jap gave a’ valué o£ 9372 Hz for 
ps, § - ce ae where J,, and J, are the coupling 


constants between the cadmium nuclei with spin numbers 


of 1/2 and H, and H,. The values for |g 





and ee. | 


AX 
shown in Table I were determined by varying them and 
calculating spectra until the best fit between 
experimental and calculated spectra was obtained (35), 
assuming the signs of the coupling constants to be 
the same. No Significant variation of these coupling 
constants with temperature was found. The stick 
spectra presented in the bottom part of Figure 3 were 
calculated from the chemical shift and coupling 
constant’s obtained above. 

The 60 MHz spectrum of the acetate methylenic 


protons of a solution of Pb (1,3-PDTA) 7 aeons? and 


pH 7.9 consists of a single resonance flanked symmet- 
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rically by two satellite peaks; the width of the central 
resonance is 1.2 Hz. The satellites arise from proton 
coupling to the ate isotope which has a nuclear 

Spin quantum number of 1/2 and a natural abundance 

of 22.6%; the magnitude of the coupling constant 

Iny is given in Table I. At both 100 and 220 MHz, 

the basic pattern of the spectrum of Pb (1,3-PDTA) “~ 
is the same as that at 60 MHz. The lack of an AB 
pattern could be due to more labile lead-nitrogen 
bonding resulting in averaging of the resonances (9) 
or to accidental equivalence of the methylenic proton 
chemical shifts (15). 

The 220 MHz spectrum of the acetate methylenic 
protons of Cd(1,4-BDTA) 77 at 5° is shown at the 
bottom of Bicune 4. The pH of this solution was 
7.6 at 25° 7 /The/spectrum consists of fhe AB part of 
an ABX pattern superimposed on an AB pattern; analysis 
of the multiplet yielded the values listed in Table I. 

The 60 MHz spectrum of the acetate protons of 
gn(1,4-BDTA) *~ at 62° consists of a simple AB pattern; 
the results of an analysis are given in Table I. 

This work was done in collaboration with Dr. G. 
Blakney (30). 
The acetate proton spectrum of Ca (EDDDA) 7” at 


220 MHz and 5° is shown at the bottom left of Figure 5. 
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Figure 4: 220 MHz pmr spectra of the acetate methylenic 


protons of an aqueous solution of 


Cd(1,4-BDTA) *~ at the indicated temperatures. 
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Figure 5: 220 MHz pmr spectra of the acetate methylenic 
é protons of an aqueous solution of 


Cd (EDDDA) 77 at the indicated temperatures. 
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It is believed to consist of the AB part of an ABX 
spectrum superimposed on an AB pattern. The two 

low intensity peaks of the AB pattern are too weak 
to be observed and hence, A6 and J could not be 


AB AB 


calculated. Values for |J and |J were estimated 


ax! Bx | 


directly from the spectrum and are listed in Table I. 
At 60 MHz and 25°, the acetate methylenic portion 
of the Zn(EDDDA) 27 spectrum consists of a simple AB 
pattern. The results’ of an analysis of this spectrum 
are shown in Table I. 
The 220 MHz spectrum of the acetate methylenic 
protons of Pb(EDDDA)* at 5° consists of the AB 
part of an ABX spectrum superimposed on an AB spectrum. 
The AB part of an ABX spectrum arises from those 
complexes containing the aon isotope... Analysis of 


this multiplet pattern yielded the results shown in 


Table I. 


AB Tibevohanweeneces for the 1,3-PDTA and 1,4-BDTA 
Complexes 

The lineshape of the multiplet pattern for the 
acetate ‘methylenic protons of Cd (1,3-PDTA) 77 is 
temperature dependent, as shown by the 60 MHz spectra 
in Figure 6 for a solution of the complex which was 


pH 8.0 at 25°. The two low intensity peaks of the 
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AB pattern are barely visible at 16° and are not 

shown in the figure. As the temperature is increased, 
the AB pattern collapses to a single, averaged resonance 
while the AB part of an ABX pattern collapses to the 

A part of an A,X pattern: The 98° spectrum clearly 
illustrates the limit of rapid interchange of the A 
and B protons of a given methylenic group on the pmr 


time scale. The coupling to the + tea and ie fe 


d 
isotopes at 98° indicates that complete dissociation 

of the complex must be occurring slowly on the pmr 

time scale, since metal-proton coupling in multi- 
dentate ligand complexes is not collapsed as long as 
chelated ligand does not exchange rapidly with free 
ligand (9). The metal-proton coupling would be observed 
even if the ligand partially unwraps from the metal 

and then re-bonds because all nuclei remain in the 

same spin states and the coupling is not relaxed. 

The upper limit of the rate of first order dissociation 


of the complex at 98° is 2r|J 87 Seca: where 


axlo* 


| J is the average value of the magnitude of the 


ax! 
cadmium-proton coupling constant in the limit of fast 
AB intefchange. The slight broadening of the 

satellite peaks at 98° suggests that at this temperature 


the rate of complete dissociation is just becoming 


appreciable on the pmr time scale (43), that is, 
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the rate must be of the order of 2m|Iny 





Average lifetimes before AB interchange were 
determined by lineshape analysis of spectra of which 
those in Figure 6 are representative. The density 
matrix method (78,79) was used to simulate spectra 
theoretically as a function of t, the average 
lifetime before interchange. Other parameters used 


J J and T 


. : H 
in the calculation were Ad.5 ( ya | ax’ Upx 21 


AB’ 
where T, is the effective spin-spin relaxation time. 
Adan at a given temperature was calculated with the 
aid of Equation I; similar chemical shift functions 
have been used previously in lineshape treatments 
(BU... £h, Constant T., calculated from an assumed 
natural linewidth of 1.0 Hz was used; the natural 
linewidth could not be measured directly because some 
exchange-broadening was observed even at the lowest 
temperatures attainable. Typical linewidths in 
similar complexes are about 1.0 Hz. 

The computer program used to simulate the spectra 
(51).was modified,so, that it.first calculated the 
lineshape of the AB spectrum and then the AB part 
of the ABX spectrum for a given set of parameters. 

It then normalized the intensities for the two 


sets of spectra so that the AB part of the ABX 


pattern represented 25.0% of the total intensity, 


. <t rts “ae ‘on ‘HK 


dtoidw to stdnege Xo eleylads wadesed yd Sohne 
Whenod off ovisedinonexdez Gee? erupit ai saode 
e1szasqe eisfumie oa? Soeu enw, (eT, aS) botten xiatam 
opexeva off ,7? Io motfonst s as viientexoelit 

beay gisiemessg, s920 _opcadoistns -sxated omisotal 
G (aH) $4 stew nokeepoodaa add ad 
omid nokieteigt nieesnide eviaoe22s of% el red otodw 





id didw bedahietés aew Siviswegmed nevig s 36 poh 
anorsonwd- Aine leolgedo. talimse, <1 noivasupa go bis 
atnemisosxt sqbdesikl at Yiguosvas® heaw seod eved 
homers Ne moss hesalucicn, oT sastande A « (82) 
forueda sc cbaetn cow an 0,129 djbiwend! Leen 
amos sausosd yitosith baxsitesen od aon BAleS nobiuentl 
Jaawol sit. .28- (ave havised> 28 oa triohsotd-spAshoxe, 
nt edsbiveanit-Lesiqy? .aitertssaa, apt Ht eTeqmea?d 
,SH 0.1 aunds-g1a 2oxekqmngg 1slimie 

sijoeaqe eid sisiumie.ot beer .msupesq zesuqmen att 
sift boteloolss devi? Ji Jed cevhettibom <av (£2) 
gus €A. orld cent Sia muadoage as ote 30 sqadeamhs 








pcm nenae ped pre i psa sr 









31 


Since this part of the pattern arises from the com- 
plexes of the isotopes of cadmium with spin numbers 
of 1/2 which are present at a combined natural 
abundance of 25.0%. The normalized spectra were 
combined to give a spectrum which was compared to 
the experimental spectrum. For a given set of 
input parameters, t was varied until the computed 
and experimental spectra matched. In the slow and 
rapid exchange regions the criterion for matching 
was the linewidth of the resonance peaks; in the 
intermediate exchange region it was the peak-to- 
valley ratio of the two central peaks. The terms 
slow, intermediate, and fast exchange will be 
explained in Chapter IV. The spectra shown in Figure 
6 are typical of the fits obtained. Average life- 
times before AB interchange determined by this 
matching procedure over a range of temperatures 
are listed in Table II. The experimental error in 
the lifetime determinations was from + 5 to £°108% 
due mainly to difficulties in matching computed and 
experimental spectra. 

A pH of 8.0 was chosen for the Cd-1,3-PDTA 
lineshape study to minimize the contributions of 
reactions involving hydrogen or hydroxide ion to the 


overall rate of AB interchange. To determine the 
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Table II 


Average Lifetimes of the Acetate Methylenic Protons 


Temp. 


16 
23 
ad 
38 
39 
45 
onl 
58 
64 


al 


(a) 


(b) 


(c) 
(da) 


of Ca(1,3-PDTA)*~ and Cd(1,4-BDTA) *~ 


Before AB Interchange 


x (sec.) 
5-arb 5c 1r4 
oot) Cd(1,3-PDTA) Cd (1,4-BDTA) 
0.10 
Om 2 
eR ote: 0.950 
0.070 
On Oe 
ie Whe pe, 
0; O42 
07 029 Bigg BY aS 
0.020 
Oe uULLE 
0.0043 


Die tee pci 

Lifetimes determined by lineshape analysis of 60 MHz 
spectra. 

Ob ii Ua teed 

Lifetimes determined by lineshape analysis of 220 


MHz spectra. 
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magnitude of any pH dependence, samples at pH's 

jen; toes ane 9, PY (at o25 ©) Were *sttdied ‘over ‘the 

same temperature range. The spectra obtained were 
identical, at the same temperatures, with those from 
which the results in Table II were obtained, 
indicating no pH dependence of the rate within this 
pH range. 

To determine if the lineshape treatment is 
dependent on the spectrometer frequency, lifetimes 
were determined by the above procedure for a D0 
solution of Cd(1,3-PDTA)* at pD 7.5 over a range 
of temperatures from both 60 and 100 MHz spectra. 
D0 was used as the solvent rather than H,O because 


2 


with HO as the solvent, the water resonance caused 


2 
locking problems on the 100 MHz spectrometer. 
The results shown in Table III indicate that the 
lifetimes are independent of the measuring frequency 
within experimental error. The results in Tables 
II and III also show that, on the average, the 
rates of AB interchange for the D,0 solution are 
0.6 times as fast as those for the HO solution. 
These results are discussed in a later section of 
this chapter. 

The 60 MHz pmr spectrum of the acetate methylenic 


protons of fx (1) 3=PDTA) 27 in aqueous solution is 
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Table, IIlt 
Comparison of Average Lifetimes of the Acetate 
Methylenic Protons of Ca (1,3-PDTA) ~~ i Before 
AB Interchange Determined from 60 MHz 


and 100 MHz Spectra 


tT Asec.4 

Temp. (°C) 60 MHz 100 MHz 
23 URiZ 

26 OL 
aL 0.10 

eh) 0.094 
45 02065 

47 OS056 
51 0.043 

56 0.034 
58 0.028 

60 0.026 


(a) D0 SERUSTONMGL TED 7S5Pat 25°. 
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shown as a function of temperature in Figure 2. 

As the temperature is increased, the AB pattern 
begins to collapse; however the extent of collapse 
and hence the rate of AB interchange is somewhat 

less “than that for Cd(1,3-PDTA) 7~, Average lifetimes 
before AB interchange were determined by matching 
computer-simulated and experimental spectra in the 
same way as for the cadmium complex as described 
above. The chemical shift difference and coupling 
constant used in the computer-simulation were found 
to be independent of temperature in the range over 
which they could be measured (25 to 62°). The life- 
times obtained are listed in Table IV, and some of 
the calculated spectra are shown in Figure 2. 

The temperature dependence of the lineshape of 
the multiplet pattern for the acetate protons of 
Cd(1,4-BDTA) 77 at 220 MHz is shown in Figure 4. 

As observed for Ca(1,3-PDTA)", the AB pattern 
collapses to a single resonance as the temperature 
is increased and the AB part of the ABX pattern is 


just beginning to become the A part of an A,X pattern 


2 
by 70°: In the intermediate interchange region, 
the lineshape is slightly asymmetric with one side 


of the multiplet sharper than the other. Average 


lifetimes before AB interchange were determined by 
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Table IV 
Average Lifetimes of the Acetate Methylenic Protons 
of Zn(1,3-PDTA)* and zn(1,4-BDTA) ~*~ 


Before AB Interchange 


t (sec.) 
ees sa 

Temps « (2C) Zn (1,3-PDTA) Zn (1,4-BDTA) 
Zo 2.9 

39 he 3 

SL Yea: Ul 

56 U. 30 

62 ge ee OF0 

Jo OL iL] 0.20 

80 a Ap i U.uss 

87 0.084 0.063 

OS 0.060 0.047 
100 0.040 0.085 


(a) Determined by lineshape analysis of 60 MHz spectra. 
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the matching procedure used for the 1,3-PDTA complexes, 
but disregarding the asymmetry. The chemical shift 
and coupling constant values used in the calculation 


are given in Table I and a T. assuming a natural 


2 
linewidth of 2.5 Hz obtained from the 220 MHz spectra 
was used. The lifetimes which were obtained are 
shown in Table II for comparison with those for 
Cd(1,3-PDTA)~. The error introduced into the 
lifetime calculations for the 1,4-BDTA complex by 
the presence of the asymmetry in the intermediate 
exchange spectra is estimated to be + 208. 

The temperature dependence of the lineshape for 
the acetate protons of gn (1,4-BDTA) 77 at 60 MHz 
is shown by the spectra in Figure 7 (30). Average 
lifetimes before AB interchange determined by matching 
computer-simulated spectra with experimental spectra 
obtained above 48° are given in Table IV for comparison 
with the results for zn (1,3-PDTA) ~*~. 

Below 48° for gn(1,4-BDTA)~, the individual 
lines of the multiplet pattern broaden at different 
rates, resulting in an asymmetrical AB pattern as 
shown in Figure 7. Each half of the spectrum is 
broadened further at 100 and 220 MHz, indicating 


that the broadening results from a second kinetic 


process. One possibility is relatively slow exchange 
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Figure 7: 60 MHz pmr spectra of the acetate 
methylenic protons of an aqueous solution 
hie Zn (1,4-BDTA) ~~ at the indicated 


temperatures. 
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of the acetate groups between the site in the plane 
formed by the metal and the two nitrogen atoms and 
the site out of the plane. Sucha process would 

6 


yield an asymmetrical pattern if J and 6. are 


AB’ A’ B 
different for the in-plane and out-of-plane acetate 
groups, as in the cobalt(III) and rhodium(IITI) 

1 i j 
conpiexes, OF? EDTA''(95 14716, 17)}" and’ Lt eA (in-plane) 


$x (out-of-plane) is different from "Bin plane 


Sei quteot-plane) It is possible that the asymmetry 
observed in the intermediate exchange spectra for 
Ca(1,4-BDTA) *~ also results from in-plane to out-of- 
plane exchange of acetate dentates. However, in this 
case, the chemical shift differences between the in- 
plane and out-of-plane sites must be very small, 


Since in the slow exchange spectrum only a slightly 


greater intensity of one of the peaks can be observed. 


femperactuLesvepenaence of the AB Multiplet Patterns 
for the EDDDA Complexes 


The multiplet pattern for the acetate methylenic 
protons of Cd (EDDDA) 77 at 220 MHz broadens asymmet- 
rically as the temperature is increased as can be 
seen by the spectra in Figure 5. At 70° it consists 
of a single resonance flanked symmetrically by the 


A part of an AX pattern due to the proton coupling 


Z 
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40 
to the ater and mireien | isotopes. The temperature 
dependence of this multiplet pattern could not be 
quantitatively studied because the chemical shift 
difference and coupling constants could not be obtained 
from the slow exchange spectra. 

The AB pattern for Zn (EDDDA) ~~ at 60 MHz broadens 
little at the highest temperatures studied, indicating 
only a slight degree of AB interchange. Lifetime 
determinations were not made because of the small 
temperature range (90 to 100°) over which broadening 
was observed. 

The multiplet pattern for Pb (EDDDA) *~ at 220 MHz 
also broadens as the temperature is increased, 
indicating AB interchange. In the intermediate 
exchange region the broadened pattern is asymmetrical. 
Severe interference from the water resonance in this 
exchange region prevented a quantitative rate study 


from being undertaken. 


Pmr Spectra of the EDTA Complexes 

The 60 MHz pmr spectrum of the acetate methylenic 
protons’ of an aqueous solution of Ca (EDTA) 77 ata 
neutral pH consists of the AB part of an ABX pattern 


superimposed on an AB pattern (9), similar to that 


‘observed for Cd(1,3-PDTA) >. In contrast to the 
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1,3-PDTA system, however, the lineshape of the 
multiplet pattern of the EDTA complex is temperature 
Hiivarzanc up to £00", 

It has previously been observed that the acetate 
methylenic proton resonances of the zinc and lead 
complexes of EDTA (9,52) at 60 MHz consist of 
Single peaks. Pmr measurements at 220 MHz (15), 
however, have revealed AB patterns for these two 


complexes. 


B. DISCUSSION 


Partial Dissociation Kinetics of the Complexes of 
1,3-PDTA and 1,4-BDTA 

The multiplets in the pmr spectra for the acetate 
methylenic protons of Cd(1,3-PDTA)*, Zn (1,3-PDTA) 7~ 
and Cd(1,4-BDTA)* at 25° and of Zn(1,4-BDTA)* at 
62° consist of single AB patterns (or the AB part 
of an ABX pattern superimposed on a single AB 
pattern). The single AB patterns indicate that the 
four acetate groups of a given complex are equivalent 
on the pmr time scale. If it is assumed that coordina- 
tion around the metal ion is octahedral and that all 
six dentates of the ligands are coordinated, two 
acetate groups are in the plane formed by the metal 


and two nitrogen atoms and two are out of the plane 
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in the instantaneous structures of the complexes. 
If a given structure is long-lived on the pmr time 
scale, the environments experienced by the acetate 
groups in the in-plane and out-of-plane sites 
must be equivalent to observe a single AB pattern. 
Poraspectra Of the ponlabile rhodium(iqI1) (14), and 
cobalt(III) complexes of EDTA (9,16,17) suggest it 
is unlikely that the in-plane and out-of-plane 
environments will be equivalent and rather that the 
Single AB patterns for the 1,3-PDTA and 1,4-BDTA 
complexes result from rapid exchange of acetate 
dentates between the two sites. Exchange could be 
by a sequence of reactions involving dissociation of 
labile metal-oxygen bonds (9) or by a twist mechanism 
of the type proposed by Bailar (18,53). 

The AB patterns reveal that, for the given solution 
conditions, each nitrogen atom of the 1,3-PDTA 
and 1,4-BDTA complexes is permanently in one of its 
two possible inversion states on the pmr time scale. 
Increasing the temperature results in inversion of 
the nitrogen atom becoming more rapid, as evidenced 
by the ‘collapse of the multiplet patterns for these 
complexes. 

A series of stepwise reactions which could lead 


to interchange of nonequivalent acetate protons in 
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complexes of EDTA and its homologs via inversion of 
a nitrogen atom is shown schematically in Figure 8. 
The ligand is shown as being hexadentate in at least 
some of the complex, but the conclusions from the 
following discussion are not dependent on this 
assumption. It is assumed, however, that a metal- 
nitrogen bond can dissociate making possible nitrogen 
inversion only after dissociation of the metal bonds 
to the two acetate groups of the given nitrogen atom. 
iba Ky (= k,/k_,) and Ky © ko/k_») are rapid 
pre-equilibria on the pmr time scale, the rate-deter- 
mining step for AB interchange via the reaction 
sequence in Figure 8 is either nitrogen inversion 

or metal-nitrogen bond dissociation. Qualitatively, 
if k_3 ates the metal-nitrogen bond dissociates 

and re-forms many times on the average before the 
noncoordinated nitrogen atom inverts, so that the 
rate-determining step is the nitrogen inversion. 

On the other hand, if Ke 7a , once the metal- 
nitrogen bond has broken the noncoordinated nitrogen 
atom inverts many times before the bond re-forms and 
the ratecdetermining step is metal-nitrogen bond 
dissociation. In the latter case, there is statis- 
tically a 50% chance that the nitrogen atom will be 


in an inverted state when it re-bonds to the metal, 
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so that the observed rate of AB interchange will be 
half the rate of dissociation to species (IV) in 
Figure 8. 

Inversion rates of nitrogen atoms in multidentate 
ligands of the type studied here are not known; 
however, they have been estimated to be 10° to 10° 
sec + for nitrogen atoms in tertiary amines of the 
type described earlier (26-28). In the absence of 
more appropriate inversion rate constants, Ks me 
nitrogen atoms in multidentate ligands of the EDTA 
type will be assumed to be of this magnitude in the 
following discussion. 

The upper limit for the rate of formation of 
the metal-nitrogen bond described by the rate constant 
k_3, is the rate of dissociation of the water molecule 
occupying the coordination site if the reaction 
occurs by the Eigen dissociative mechanism (31-33,70). 
In terms of this mechanism, a water molecule will 
first be released from a site cis to the site of 
coordination to the other nitrogen atom, and the 
noncoordinated nitrogen atom will then bond to the 
metal abe at the vacated coordination site. The 
rates of water release from totally aquated cadmium 


8 ¥, 


Bnd Pixc Sr err ytg 2 ana bem §4e, 687+ respectively 


(54). The rates of water release from these metal ions 
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will be considered as the upper limits for k_3 
in the following discussion. 

For k_3 to equal the rate of water release, the 
rate at which the noncoordinated nitrogen atom in 
species (IV) of Figure 8 coordinates to an unoccupied 
Site has to be significantly larger than the rate at 
which a water molecule reoccupies the site. This 
requires that the nitrogen atom be properly oriented 
for bonding at the time the water molecule dissociates. 
Previous investigations of a number of fast complexa- 
tion reactions of multidentate ligands would suggest 
that the rate at which the nitrogen atom will coordinate 
to the unoccupied site will be less when the chelate 
ring formed is six-membered than when it is five- 
membered. For example, Kustin and coworkers observed 
that the formation rate constants of the g-alanine 
complex of cobalt(II) is almost one order of magnitude 
less than that for the a-alanine complex (47). The 
difference was attributed to the difficulty in forming 
the six-membered chelate ring. Similar differences 
were observed for the nickel and cobalt(II) complexes 
of a- ahd g-aminobutyric acids (48). These results 
suggest that, for Cd(1,3-PDTA) “7 and zn(1,3-PDTA)~, 


k may be somewhat less than the rates of water 


=o 


release from species (IV) in Figure 8 due to the 
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formation of a six-membered chelate ring. If so, k; 
would be somewhat greater than k_3 for these complexes. 
On this basis, it is proposed that the AB patterns 
observed for these complexes are indicative of 
relatively inert metal-nitrogen bonding on the pmr 
time scale, and that the inverses of the lifetimes 
before AB interchange, which are listed in Tables II 
and IV, are half the first order rate constants for 
partial dissociation of the complexes to species (IV) 
at the indicated temperatures. If this proposal is 
correct, the lifetimes before AB interchange indicate 
that the relative order of lability of the metal- 
nitrogen bonds in the two metal complexes of 1,3-PDTA 
tS CacmiUuM 2 FLAC. 

The rate at which the metal-nitrogen bond in 
the 1,4-BDTA complexes dissociates should be the 
same or slightly faster than in the 1,3-PDTA complexes 
since it may be easier to open the larger seven- 
membered chelate ring (47,48), while the rate at 
which the bond re-forms is predicted to be less 
since a seven-membered ring is formed (1,58). There- 
dapac yy te k; is somewhat greater than k_, for the 
1,3-PDTA complexes, the rates of interchange in the 
1,4-BDTA complexes of cadmium and zinc should be 


the.same or slightly faster than in the 1,3-PDTA 
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complexes. The results in Tables II and IV indicate 
this to be the case. Furthermore, the lifetimes 
before AB interchange indicate that the order of 
lability of the metal-nitrogen bonds for the complexes 
of 1,4-BDTA is cadmium > Zinc. 

To determine if the rate constants calculated 
from the lifetimes in Tables II and IV are reasonable 
for the reaction represented in Figure 8, comparison 
can be made with the rate constant Kar for dissociation 
of one N-methyliminodiacetic acid (MIDA) ligand 


from the complex Zn (MIDA) 5”, 


CH,COoH 
Coa 


CHACOln 


MIDA 


k 


= d nie 
zn (MIDA) 5 a GOLA) MIDA 


which has previously been determined by pmr line- 
broadening measurements (30). The value for Ka is 
AGez*S ae cot at 25°. Taking the inverse of the life- 
times before AB interchange to be half the first order 
rate constants for partial dissociation to:species (IV) 


in Figure 8, the partial dissociation rate constant 
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for gn(1,3-PDTA) 7~ is 40969 Bee A see ta 25, 
The difference between k, for zn (MIDA) 5” and this value 
is of the magnitude expected if the energy barrier 
associated with the necessary twisting of the ali- 
phatic chain in Zn(1,3-PDTA)2- is accounted for (1). 
The results in Tables I and II which indicate 
that the rate of AB interchange in Cd(1,3-PpTA) 77 
is different in HO and D9, can be accounted for if 
the fraction of the complex that exists in the various 
partially coordinated forms, shown in Figure 8, is 
different in these two solvents. This would be the 
case if the rate at which HO and D0 dissociate 


from the metal complex is different. For example, 


Ky ( = kj /k_}) would bejdifferent in HO and DO ix 
the mate ofidissociation of HO from species (IT) 
is/diftferent from the rate dissociation ‘of DO. 


Temperature Dependence of the AB Multiplet Patterns 
of the EDDDA Complexes 

The multiplet patterns observed for the acetate 
methylenic protons of the cadmium, zinc and lead 
complexes of EDDDA indicate an inertness on the pmr 
time scale of the metal-nitrogen bonds. The collapse 
of these patterns with increasing temperature indicates 


that partial dissociation of these complexes is 
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occurring on the pmr time scale. The different 
degrees of collapse of the multiplet patterns for 
these complexes suggest that the relative order of 
lability of the metal-nitrogen bond is cadmium > lead > 
Zinc 

Assuming that coordination in the EDDDA complexes 
is octahedral and that all six dentates are coordinated 
to the metal ion, three instantaneous structures are 
possible depending on whether the acetate or 


propionate groups are in-plane or out-of-plane (59). 


QO N 0 ~ 0 N 
0 N 0 hg ©) : 
pet: ys re 


(T) (I) (IIT) 


The single AB pattern observed in the pmr spectrum 
for Zn (EDDDA) ~~ suggests rapid exchange of the 
carboxylate dentates, but the values for Jap and Adan 
indicate that one of the above structures is pre- 
dominant. The magnitude of the geminal coupling 


constant Jap is approximately the same for all the 
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51 
complexes listed in Table I except that for Zn (EDDDA) ~~. 
Sudmeier and coworkers (17) have reported the magnitude 
of the geminal coupling constants of the methylenic 
protons of the in-plane and out-of-plane acetate groups 
in the inert cobalt(III) complexes of EDTA type 
ligands’ to be 16 +°0:5 Hz and 18 t 0.5 Hz respectively. 
Since geminal coupling is sensitive to the electro- 
negativity of substituents, the magnitudes of in-plane 
and out-of-plane coupling constants of aminocarboxylic 
acid ligands coordinated to divalent metal ions may 
be different. Assuming the relative magnitudes to 
be in the same direction as those for the cobalt 
complexes, the results in Table I suggest that 
gn (EDDDA) “> exists predominantly as structure (IIT) 
shown above. Support for this structure also comes 
from the observation of Legg and coworkers (60) that 
Adan for out-of-plane acetate groups of cobalt(III)- 
aminocarboxylic acid complexes is larger than for 
in-plane groups, and the results of Byers and 
Douglas (59) which show the structure of one isomer 
of Co(EDDDA) to be similar to the one proposed 
here for Zn (EDDDA) ~. 

If the preceding conclusions are correct, the geminal 


2=- 
coupling constant reported for the:AB pattern of Pb(EDDDA) 


as the temperature is increased might then be explained 
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by a temperature dependence of the populations of 
the above three structures. The value for Jap could 
not be determined from the multiplet pattern for 
Ca(EDDDA)*", but the asymmetric broadening of this 
pattern might also be explained by a temperature 
dependence of the populations of the given structures. 
The lack of broadening in the multiplet for 
Ca (EDTA) 77 up to 100° indicates that the rate of AB 
interchange is significantly less than the rate in 
Cd(1,3-PDTA)*>. This difference could result from 
more inert cadmium-nitrogen bonding in Ca (EDTA) 77, 
as predicted by the observation that it is easier 
for a six-membered chelate ring to open than for a 
five-membered ring (1,58), or it may indicate that 
k_3 >> k;. Bond re-formation, the reaction step 
described by the rate constant k_3) is predicted 
to be faster for Cd (EDTA) 77 than for the 1,3-PDTA 
complex since a five-membered chelate ring is being 
formed (47,48). Considering the EDDDA results, 
the rate of partial dissociation to yield species 
(IV) in Figure 8 is predicted to be larger for 
Cd (EDDDA) 77 than for the EDTA complex on the basis 
that intermediates in the stepwise sequence will 


be present at higher concentrations while the rate 


of cadmium-nitrogen bond re-formation is assumed to 
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be the same. Thus, if the AB pattern of Cd(EDTA) 
is due to the cadmium-nitrogen bond re-forming before 
inversion occurs, the AB pattern of Cd (EDDDA) 77 
should be temperature invariant. Rather the multi- 
plet for the acetate methylenic protons of the EDDDA 
complex is completely collapsed to a single resonance 


frankea "by the A ™partvor an Alyx" pattern at’'70", 


2 
indicating rapid AB interchange on the pmr time scale. 
These considerations provide support for the proposal 
of Day and Reilley (9) that the multiplet pattern 


for Ca (EDTA) 77 results from inert cadmium-nitrogen 


bonding. 
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CHAPTER IIT 
LIGAND EXCHANGE KINETICS OF METAL COMPLEXES 


OF 1,3-PDTA AND EDDDA 


In this chapter, the results of a pmr line- 
broadening study of the ligand exchange kinetics 
of the 1,3-PDTA and EDDDA complexes of cadmium, zinc 
and lead are presented. These systems were investi- 
gated with the objective of elucidating the mechanisms 
of ligand exchange reactions of cadmium, zinc and 
lead complexes of aminocarboxylic acid ligands with 
emphasis on the formation reactions of the mono- 
protonated forms of the ligands. These two particular 
ligands were chosen in order to observe the effects 
of slight structural changes on the rates of reactions 


of aminocarboxylic acids of,the EDTA type. 


A. RESULTS 
pH Dependence of the Chemical Shifts 

The pmr spectrum of 1,3-PDTA consists of a sharp 
singlet for the acetate methylenic protons, a broad 
triplet;like pattern for the protons of the two 
nitrogen-bonded methylenic groups of the propylene 
part of the ligand and a broad quintet-like pattern 
for*theicentral methylenic protons over the pH range 


2 to 13. The lack of any broadening of the exchange- 
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averaged acetate methylenic resonance indicates that 
exchange of the acidic carboxylate and nitrogen protons 
is fast on the pmr time scale. The resonances due 
to the propylene protons, which form an A,B spin 
system, were not used in the ligand exchange studies 
because of their complexity. 

The chemical shifts of the 1,3-PDTA resonances 
are pH dependent as illustrated by curve A for the 
acetate methylenic protons in Figure 9. By analogy 
with EDTA, the two nitrogen atoms of the tetraanion 
of 1,3-PDTA are protonated first, followed by 
protonation of the carboxylate groups (61-63). From 
the chemical shift vs. pH data, values were determined 
for the third and fourth ionization constants of 
1,3-PDTA, defined by Equations II and III, by methods 


described previously (44,67) and are given in 


a [HL] 
K,=— y) =7 
3 = 
[H5L ] 
ay uead 
i See Bit 
[HL J 


Table V. These constants are mixed mode constants 
involving hydrogen ion activity and 1,3-PDTA con- 


centrations. 
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CHEMICAL SHIFT, ppm vs TMA 


Figure 9: 





pH dependence of the chemical shift of the 
acetate methylenic protons of 1,3-PDTA. 


Curve A, 0.20 Mol, 3-PoTA; curve B,0.20 M 


Cd(1,3-PDTA) 7 s Curves, 0.20 M ZniL, o~ 
2=- 


PDTA)*"; curve D, 0.20 M Pb(1,3~PDTA) 


56 


# Io gRitle Ismimodo eds 30 sostebasaeb Bq 





pK) 
PK, 
PK, 


PK, 


(a) 
(b) 


(i) 


Table V 


Acid Ionization Constants of 1,3-PDTA 
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Over the pH range 3 to 13, the pmr spectrum of 
EDDDA consists of a sharp singlet for the acetate 
methylenic protons, a sharp singlet for the ethylenic 
protons between the two nitrogen atoms and a complex 
multiplet pattern for the carbon-bonded protons of 
the two propionate groups. The spectrum at vH 10 
has been described previously (59). On the basis 
of chemical shift data for other structurally- 
related aminocarboxylic acids (62,63,68), the down- 
field singlet was assigned to the acetate methylenic 
protons. The chemical shifts of the EDDDA resonances 
are pH dependent as illustrated by curve A for the 
acetate methylenic protons in Figure 10. The third 
and fourth ionization constants were determined from 
these data; the results are given in Table V. 

The spectra of the acetate methylenic protons of 
the cadmium, zinc and lead complexes of 1,3-PDTA 
and EDDDA are revorted in Chapter II. In most cases 
they consist of complex multiplet patterns due to 
inert metal-nitrogen bonding. To a first approximation, 
the resonances of the remaining carbon-bonded protons 
of tess two complexed ligands appear the same as 
those of the uncomolexed ligands, but have different 
chemical shifts. 


The chemical shifts of the acetate methylenic 
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Figure 10: pH dependence of the chemical shift of 
the acetate methylenic protons of EDDDA. 


Curve A, 0.20 M EDDDA; curve B, 0.20 M 
Ca(EDDDA)*-; curve C, 0.20 M Zn(EDDDA) 7"; 


curve D, 0.20 M Pb (EDDDA) ~~. 
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protons of the cadmium, zinc and lead complexes of 
1,3-PDTA and of EDDDA are pH dependent as shown by 
curvyes«<By ,-Ceand*DeinSPigures” 9° and- 10.°*'These data 
were obtained from solutions made up from equimolar 
quantities of both metal nitrate and ligand. For 
those complexes in which an AB pattern was observed, 
the chemical shift values used in these curves were 

the chemical shifts at the centres of the patterns. 

For all the metal complexes studied, the chemical 
shifts remain constant above a certain pH, corres- 
ponding to the chemical shifts of the 1:1 complexes. 
Below this pH, the resonances move towards those of 

the free ligand and are exchange-broadened. Calculations 
of the species present in solution for the 1,3-PDTA 
systems using the literature ionization constants in 
Table V and the literature formation constants 

(64-66) indicate that dissociation of the 1:1 complexes 
and formation of protonated complexes is occurring 

at those pH's when the resonances are shifting and 

ind broadened. For example, calculations indicate 

that ina solution with a cadmium to 1,3-PDTA ratio 

of snefedi sedeudeten of the 1:1 complex and formation 
of a protonated complex begin to take place at pH 

less than 5.5. It is also below this pH that broadening 


and shifting of the resonance begins tc take place 
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as can be seen from Figure 9. For the EDDDA systems, 
no formation constants are available for the protonated 
complexes, but such species are indicated by the 
exchange-broadening and the chemical shift data shown 
in Figure 10. 

An attempt was made to calculate formation constants 
for the nonprotonated and the protonated complexes 
from the chemical shift data presented in Figures 9 
and 10 using the method described previously (44). 
These formation constants would have been useful in 
subsequent calculations since they would have been 
obtained under ionic strength conditions similar to 
those used in all experiments reported in this work. 
This attempt, however, yielded inconsistent results, 
possibly because of the uncertainties in the chemical 
shift values of the observed resonances and those of 
the protonated complexes. The observed resonances 
are broad and hence it was difficult to measure 
their chemical shift positions. The chemical shift 
values of the protonated complexes could not be 
measured directly, but had to be approximated in the 
attempted calculation of formation constants. 
Therefore, in the subsequent rate constant deter- 
minations it was necessary to use literature values 


of the formation constants (64,65). However, the 
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error involved using these literature values is 
within the experimental error of the lifetime measure- 


ments. 


Ligand Exchange Kinetics of the 1,3-PDTA Complexes 


Exchange-broadening was observed in the pmr 
spectra of each of the 1,3-PDTA complexes under 
conditions where some noncomplexed ligand was present. 
From the broadening for a given set of experimental 
conditions, the mean lifetime t between events which 
lead to exchange of the ligand from the free form 
to the complexed form or vice versa may be determined. 
The following reactions are considered to be pathways 


by which exchange can occur: 
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63 
methylenic protons of Ca (1,3-PDTA) *~ and zn (1,3-PDTA) >~ 
at intermediate pH values in the absence of non- 
complexed ligand, which were discussed in detail 
in Chapter II, indicate that for these conditions 
the mean lifetimes before AB interchange are greater 
than 0.1 sec. This value was obtained by computer- 
Simulating spectra by the method described in Chapter 
II. Since AB interchange may occur even when the 
complex is ‘partially bonded, the above lifetime is 
a lower limit for the mean lifetime of the complex 
before it completely dissociates via the reaction 
represented by Equation IV and indicates that the 
rate of this reaction is too slow to be measured 
by the pmr line-broadening method. Similarly the 
metal-proton coupling observed for the cadmium and 
lead complexes in the absence of noncomplexed ligand 
indicates that for these complexes the reaction 
represented by Equation IV contributes a negligible 
amount to the exchange-broadening (43). Neglecting 
this reaction, the rate equation for exchange of 


ligand from the free form to the complexed form is 


+ 


~d{Le] 2 3 


= k_5(M?*) [HL77) + k,(L 1 [ML ~] + k, (HL MiMi 


dt 


2- a 
r k, [HoL J] (ML™ J 1.4 






Ciera mt bescuerb ovew cole s onspit bexetamoy 
gaokjibmoo eeott zo? Sent eteskbak y TT. wstqedd nt | 
1etse1%% o76 spnstinaodal dé orcTed sombbeRet dager edt 
_rodiqmos yt beabesdo eow opkae eee yen Due edd 
<siond? ‘i fediteesh bodtom ale of Bassge tpaeee tome 
oi? notw asve asso yer opandonetal BA soaks> oer 





si omttetil svedse oe \bebned yilsleaeqe? xpigeese 


xoltqnos edd Yo omtjetit mgem ent 202 simkt teweks 


"wotjoser afd siv esteksozerh visseiqnes 2% eseied 


sii tndt sotsoibat bom VT oideupd yd  Sesasesaged 
bothesom od OF wola cos ak nolsseetp aidd Io tex 
sit yfxslinis s»borttem orinisbsordsentl amp eatygs 
bas muimbss edd 16% Bevesedo prifaes nosorg-Iston 
boepit Pexedgmosonon io spieads sit gi aueeianes Baek 
noijnss1 sd eexsfdmon saane ac? Sede sepa kiat 
alfdipitpen s estddixzdqeo VI acitdenpy yd begneeedees 
enidoalpen .(£&) poineshsord-stesiose eff GF Saas 
to senetioxs vol gaiasaps 6252 Se algouas eid 
ei mot baxelgqnoo oils oJ wxo2 ostt edd moxt heepet 
7 ; ie ' F 








64 


where [L-] represents the sum of the concentrations 
of the free ligand in the various protonated forms. 
The reciprocal of the mean lifetime of the free form, 


1/ty , is related to the exchange rate, -d{L,]/dt, by 
£ 





which yields 





Loe 2+ 2- 2- 2- 
= = k_5%31M ] + ka, [ML ] + k 03 [ML } + k,a. [ML ] 
£ 
MAE 
ae ce dense [HUD whe Goel SA eee /toeeane 
2 2 cseatcinan bam aeiagaateg 
a, = [Lenk /dler. Expressing the concentrations in 


terms of fractional concentrations and the term 
for the contribution from the reaction represented 


byihguation Vani termsdok [ML] by the relation 


Z 


k_5(M**] = k,[H"] (ML* ]/a,[L,] yields 


a + 
) = k, [H ) + (k + kia, + Keay) [Le] Xi 


Bn (5 i 304 a5 


Le ML 


where Py and PML are the fractional concentrations 
fi 


of ligand in the free and complexed forms (PL + 
£ 


PuL = 1). The individual rate constants were evaluated 


from the dependence of tT, on (HJ and [L,]. 
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A typical example of the spectra obtained as a 
function of pH for a solution 0.10 M in both Cd(1,3- 
PDTA) 77 and free 1,3-PDTA is given in Figure ll. 

At pH 6.42 separate acetate methylenic proton 
resonances are observed for both the free and the 
complexed forms of 1,3-PDTA; only the two strong 
central peaks of the AB pattern for the complexed 
form can be seen. An increase in the pH results in 
broadening of the resonances indicating faster 
exchange. Since the linewidth of the free ligand 
resonance for most of the systems in this study 
could easily be measured over a wide range of solution 
Conditions, ine was calculated using Equation 
XVIII, the slow exchange expression, which will be 
discussed in more detail in Chapter IV. The errors 
involved in determining the values of 1/ty arise 


i 
mainly from the linewidth measurements and are 


estimated to be between + 5% and + 158. 

Lifetime data are given in Table VI for the 
Cd-1,3-PDTA system ae experiments in which the pH 
was held constant at six different values and the 
concentfation of free ligand was varied at each 
value, »A typical plot mwsing the data at.pH 6.51 45 
shown in Figure 12. According to Equation XII, the 


slope of a plot of PL PAG: vs. [Le] at constant 
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pH 7.14 
PH 6.42 
~N-CH,-C 0; 2N-CH,~CO, 
FREE COMPLEXED 
Figure 11: 60 MHz pmr spectra showing the acetate 


methylenic resonances for an aqueous 
solution containing both free 1,3~-PDTA 
and Cd(1,3-PDTA)2- at the indicated pH 
values. The solution was 0.1 M in both 
free 1,3-PDTA and Cd(1,3-PDTA)2-. 
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Table VI 


Kinetic Data for the Cd-1,3-PDTA System at Constant 


PH Values as a Function of Free Ligand Concentration? "3 
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(Table VI, continued) 
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PH is equal to K3%, + Ko + kro, and the intercept 
is equal to ke [oe ls At the pH values listed in 
Table VI, “, ranges from factors of 10° Los 10° 
greater than ay while k3 has been found to be only 
10 to 20 times larger than kK, in similar reactions 
of EDTA complexes (39,40), suggesting that the term 
ko, contributes a negligible amount to the slopes 


of the ipa a 2h vs. [Lp] plots. Assuming the 


+ ko 


slope m to be equal to k 3 5%" 


a Equation XIII 


4 


can be derived. 


m/a, = k, (43/45) + k BCT ar 


2 5 


Rate constants ky and ke were evaluated from the 
slope and intercept of a plot of m/ 4 vs. 15/A5; 
the values for O5 and a, were calculated from the 
literature dissociation constants listed in Table V. 
A linear least squares analysis was used for this 
and all subsequent determinations of slopes and 
intercepts. The rate constants and their standard 


deviations are given in Table VII. The uncertainty 


in the ri te |e of the plots of PL Bite Pur) vs. 
E 


[L, ] ar which the data are given in fie VI, was 
too large to allow evaluation of k,- 
Lifetime data are given in Table VIII as a 


function of pH for a solution 0.10 M in both 
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Table VIII 


Kinetic Data as a Function of pH for the 


pH 


(a) 
(b) 


eet 


Cd-1,3-PDTA system?!?'° 


Wise2 (Hz) Man. (sec) 
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Cd(1,3-PDTA) *~ and free 1,3-PDTA. As can be seen 
from these data and from the spectra in Figure ll, 
the rate of ligand exchange decreases as the pH 
is increased from 4.8, goes through a minimum at 
pH 6.4 and then increases as the pH is increased 
further. Above pH 8.5, the exchange rate is so 
fast that the lineshapes are of the intermediate 
to fast types. Since exchange occurs between a 
single line, an AB pattern, and the AB part of an 
ABX pattern, the theory is quite complicated and 
no attempt was made to evaluate lifetimes from data 
at pH's greater than 8.5. 

Rate constant ky was obtained from the data in 
Table VIII at pH less than 6.4. At these pH's, 


ay is negligible relative to a. and a and the 


2 3 
contribution to the total exchange by the reaction 
represented by Equation VI may be neglected. Since 
Pp equals P 

L 
oi 
to 


ML’ rearrangement of Equation XII leads 


z 





= + 
- (kK, %3 + koao) [Le] = k, [H ] Ly. 


T 


5 a 


The value given in Table VII for ko was obtained from 


the slope of a plot of / bug - (K,a3 + keay) [Le] vs. 


(H’] using the data in Table VIII and the values 
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determined above for Ky and Kee This plot is shown 


in Figure 13. Rate constant k_5 was obtained from 


k, using the relation k_, = k,K,K, where K, is the 


2 eae as f£ 
formation constant of Cd(1,3-PDTA) 77 (64-66). 

Rate constants Ky and kK, were also determined 
from the data in Table VIII in the pH range 6.4 to 
8.2. In this pH range, calculations using the value 


of k, show the amount of ligand exchange via the 


Z 
reaction represented by Equation V to be negligible. 
The reaction represented by Equation VI has also been 


neglected for reasons given above, leading to 





iu O. 
= Sam mee peas As + k XV 
ty Pel %5 Ae ae s 
. 
by rearrangement of Equation XII. The values obtained 


£Gr ky and ke from the slope and intercept of a plot 
of Ts Bh} vs. 13/0. are given in Table VII. 

Above pH 8.5, the rate of ligand exchange 
increases, as evidenced by the intermediate to fast 
exchange spectra, presumably due to increased con- 
tributions from the reactions represented by 
Equations VE ®and “Vr 3 

tifetimd data are given in Table IX as a function 
of pH for a solution 0.10 M in both zn (1,3-PDTA) ~~ 


and free 1,3-PDTA. Slow exchange spectra were obtained 


throughout the whole pH range studied and Equation 
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Table IX 


Kinetic Data as a Function of pH for the 


pH 
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XVIII was used to calculate the lifetime Th from 
the acetate methylenic proton resonance of the free 
1,3-PDTA. As in the Cd-1,3-PDTA system, the rate 
constants ky and Ke were obtained from the slope and 
intercept of a plot of ao Barer saet vs. a3/ 45 using 
the kinetic data at pH's greater than 7. The value 
for ky was obtained from the slope of a plot of 
es 
Le 
data at pH less than 6 and the values of Ky and ke 


determined above. All rate constants so obtained 


(k,03 + Keay) [Le] vs. [H’] using the lifetime 


are listed in Table VII. 

Intermediate to fast exchange spectra for the 
acetate methylenic protons are exhibited by the 
Pb-1,3-PDTA system; examples of such spectra are 
shown on the left hand side of Figure 14 for a 
0.10 M solution of Pb(1,3-PDTA)*~ at pH 7.00 which 
contained the indicated concentrations of free 1,3- 
PDTA. Because of the complexity of this exchange 
case, kinetic data were obtained from the experi- 
mental spectra by comparison with computer-simulated 
spectra. The simulation was accomplished with a four- 
site exchange program (43) formulated using modified 
Bloch equations. The four sites are free 1,3-PDTA 
and 1,3-PDTA complexed to lead in three different 


207 


nuclear spin states, +1/2 and -1/2 for Pb and 
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zero for all the other lead isotopes. The procedure 
was to simulate a series of spectra as a function of 
the mean lifetime of the Pb-1,3-PDTA complex before 


ligand exchange t using the following input para- 


ML’ 
meters: populations, chemical shifts and natural 
linewidths at each of the four sites. The individual 
lifetimes for given solution conditions were then 
obtained by comparing the linewidths of the resonance 
peaks and by superimposing the experimental and 
calculated spectra. 

The lifetime data presented in Table X were ob- 
tained from experiments in which the pH was held 
constant and the free ligand concentration was varied. 
The matched experimental and computer-simulated 
spectra shown in Figure 14 are typical of the fits 
obtained in such experiments. Applying Equation XII 
to the above lifetime data, the slopes of the plots 


OLeE vss [Le] yielded the values for k3, 


(a P 
Ly Le ML 
and (ky -listedsin table Vil... The intercepts of such 


4 5 
plots gave an approximate value of 2.2 x 107 Masse 


k 


for Ky: 


A riiore precise value for ky was obtained using the 
data given in Table XI which gives the pH dependence 
of 1/1, for a solution 0.10 M in both Pb(1,3-PpTA) 


f£ 
and free 1,3-PDTA. Equation XIV was applicable in 
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Table X 


Kinetic Data for the Pb-1,3-PDTA System at 


Constant pH Values as a Function of Free 


Ligand Concentration S 
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Table XI 


Kinetic Data as a Function of pH for the 


pH 


oR oo § 


(a) 
(b) 
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Pb-1,3-PDTA System 2'P 
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this case, so that the slope of a plot of 
1/t, ~- (kya, + k,a,)[L,] vs. [H'] yielded the value 
for k, shown in Table VII. Rate constant k_5 was 
obtained from the relation k_» = kK,K.. 
Ligand Exchange Kinetics of the EDDDA Complexes 

The reactions represented by Equations IV to VIII 
were considered to be the predominant pathways by 
which EDDDA exchanges between the free and complexed 
forms fOr ,thessolution yconditions ,»used in this work. 
As in the 1,3-PDTA systems, the AB multiplet patterns 
observed for the acetate methylenic protons of the 
EDDDA complexes at intermediate pH values in the 
absence of noncomplexed ligand indicate that the mean 
lifetime of the complexes before first order dis- 
sociation by the reaction represented by Equation IV 
is greater than 0.1 sec. Thus, this reaction can be 
neglected in the kinetic analysis of the data. Free 
ligand lifetimes were determined from the linewidths 
of the exchanged-broadened acetate resonance of the 
free ligand using Equation XVIII. 

Lifetime data for the Cd-EDDDA system for experi- 
ments in which the pH was held constant and the free 
ligand concentration was varied are given in Table 


XII. Using Equation XII, the slopes obtained from 
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Table XII 
Kinetic Data for the Cd-EDDDA System at 
Constant pH Values as a Function of 


Free Ligand Concentration PrP 
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plots of Py TA 3 yievs'. [Le] yielded the value for 
£ 


L, ML 
k3 listed in Table VII and values of zero within 
experimental error for Ky and key indicating that the 
latter rate constants are too small to be measured 

by the pmr method. From the lifetime data as a 

function of pH for a soluticn 0.10 M in both Cd (EDDDA) 77 
and free EDDDA shown in Table XIII, the value for k, 
was determined. Over the pH range studied, Equation 


XII could be reduced to 





+- 
= k, [H ] XVI 


The terms in ky and Ke were neglected because it was 


shown above that their exchange contributions are too 
small to be detected by pmr line-broadening methods. 


Knowing the value of k calculations indicated that 


ce 
the term containing this rate constant could also be 


neglected. For the given solution conditions 


Be ne 


[H ] yielded the value for k 


and the slope of the plot of 1/ty, vs. 
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> shown in Table VIL. 


Lifetime data for the Pb-EDDDA system at constant 
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pH'varying the free ligand concentration are presented 
in Table XIV. Using Equation XII the slopes obtained 
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Table. XPLs 


Kinetic Data as a Function of pH for the 
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Table XIV 
Kinetic Data for the Pb-EDDDA System 


at Constant pH Values as a Function of 


Free Ligand Concentration a 
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experimental error for ky and ke. Data as a function 


of pH for this system are given in Table XV. In 


the determination of Ko from these data, Equation XII 


was reduced to 





: = ate 
: = K3%), [Le] = k, [H ] XVII 


In this case the term containing k, was too large 


3 


to be neglected over the pH range shown in Table XV. 


+ 
Therefore, a plot of 1/t, ~ k30, [Le] Veen | Sa 


£ 


yielded the value of k. given in Table VII. 


2 
For the Zn-EDDDA system it was found that the 

rates of the reactions represented by Equations VI 

to VIII are too slow to be measured by pmr line- 


broadening methods. The value for k however, was 


ha 
determined from the slope of a plot of 1/ty vs. 


E 
[H"] using the lifetime data in Table XVI. 


BB, DISCUSSION 
Mechanism of Formation of the M(1,3-PDTA) 77 and 
M(EDDDA) *~ Complexes 

In the present investigation, the rate constants 
for the complexation of cadmium and zinc by mono- 
protonated 1,3-PDTA (k_» in Table VII) were found 


to be essentially the same as those with mono- 
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Table XV 


Kinetic Data as a Function of pH 


for the Pb-EDDDA System 2/P/ 
W! ,. (Hz) 1/t, (sec™*) 
1/2 Ly 
3.08 6.15 
2.50 4.34 
Bee 3.62 
F256 3.27 


2586. 
0.10 M Pb(EDDDA)*~ and 0.10 M free EDDDA. 
Wi /2 of free ligand resonance for no ligand 


exchange is 1.03 Hz. 
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Table XVI 


Kinetic Data as a Function of pH 


for the Zn-EDDDA System cena 
' a, 
pH Wa 2s8z) ene eee ) 
Se es: aks hey 6.28 
= eg cy 2.46 4.62 
2 1.88 2.80 
4.06 A liye be) 1 hg oh 
4.26 a eat Rate 
4.46 ea 0.54 


CA li prea re 
(b) 0.10 M Zn(EDDDA) ~~ and 0.10 M free EDDDA. 


(c) W of free ligand resonance for no ligand 


jie 


exchange is 1.03 Hz. 
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protonated EDTA (15,24,34). This finding contrasts 

with that of Pearson and de Wit (69) who determined 

that the rate of the reaction of zinc with monoprotonated 
EDTA was less than that with monoprotonated 1,3-PDTA 
using stopped-flow kinetic measurements. Their 

results can only be considered approximate, however, 

Since the observed rates were just at the detection 

limits for their stopped-flow apparatus. The rate con- 
stants for complexation of cadmium, zinc and lead by mono- 
protonated EDDDA are from 2.5 to 10 times smaller than the 
corresponding ones with monoprotonated 1,3-PDTA. 

In the dissociative mechanism (1,70) of complex 
formation, the reaction takes place in a stepwise 
fashion, with metal ion and ligand first diffusing 
together to form an outer-sphere complex. A water 
molecule then dissociates from the aquated metal ion 
and a donor group of the multidentate ligand coordinates 
at the vacated site. When a second water molecule 
dissociates” at "a position Cis’ to the metal—Ligand 
bond the second metal-ligand bond forms. The reaction 
continues with a separate step for each of the bonds. 
Depending’ on the structure of the ligand, the rate of 
formation of a specific metal-ligand bond will 
govern the complexation rate. For the majority of 


complexes whose formation kinetics have been studied, 
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formation of the first bond is rate-determining, its 
rate being governed by the rate of water dissociation 
from the aquated metal ion in the outer-sphere complex 
CdD Ive 

A possible sequence for the stepwise reaction of 
aquated metal ions with monoprotonated hexadentate 
ligands of the EDTA type is shown in Figure l. 
Reaction proceeds by dissociation of water molecules 
followed by metal-ligand bond formation to yield 
an intermediate in which the three dentates of one 
end of the ligand are metal-coordinated while 
coordination at the other end is blocked by the proton 
on the nitrogen. For further coordination to occur, 
the second metal-nitrogen bond has to form before 
the intermediate dissociates to reactants. Presumably 
the second metal-nitrogen bond can form only after 
the proton has been transferred to a carboxylate oxygen 
er tora solvent molecule. 

The similarity of the rate constants for the 
reaction of monoprotonated EDTA and 1,3-PDTA with 
aquated cadmium and zinc (k_») suggests that, once 
the nitrogen-protonated intermediate has formed with 
these ligands, coordination of the remaining dentates 
is not rate-determining. If it were rate-determining, 


it is likely the overall complexation rate constants 
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for monoprotonated 1,3-PDTA would be less than those 
for monoprotonated EDTA due to closure of the six- 
membered ring which is necessary in the former 
complexes. Kustin and coworkers (47,48) have found 
that when ring closure is the rate-determining step, 
the rates of certain complex formation reactions are 
slower when six-membered rings are involved than when 
five-membered rings must be closed. The present 
results are consistent with the rate-determining step 
being formation of one of the first three metal- 
ligand bonds, its rate being governed by the rate of 
water loss from the metal ion. 

The smaller rate constants for reaction with mono- 
protonated EDDDA (k_>) suggest that the rate-deter- 
mining step is different in these reactions. It is 
possible that, when the first ligand dentate to 
bond is a propionate carboxylate oxygen, the rate 
of formation of the metal-nitrogen bond to close 
the six-membered chelate ring is less than when the 
first ligand dentate to bond is an acetate carboxylate 
oxygen. If this were the only factor responsible 
for the*smaller rate constants, the rates of complex- 
ation with monoprotonated EDDDA should be at least 
half those with monoprotonated 1,3-PDTA. 


A larger decrease, the magnitude of which would 
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be metal ion dependent, might result if the rates 

of dissociation of the nitrogen-protonated inter- 
mediates in the EDDDA reactions are sufficiently 
large that a significant amount of intermediate 
dissociates to reactants prior to proton transfer 
from the nitrogen atom. Kinetic data are not available 
for model tridentate complexes from which the rates 
of dissociation of the nitrogen-protonated, tri- 
dentate intermediates of 1,3-PDTA and EDDDA can be 
obtained. However, assuming the rate constants for 
formation of the nitrogen-protonated intermediates 

to be equal to the rate constants for the reaction 
with monoprotonated 1,3-PDTA (k_,) and the formation 
constants for the intermediates to be approximately 
those of complexes of model tridentate ligands, order 
of magnitude rates of intermediate dissociation can 
be predicted. The dissociation rate constant is 

¢! where Ke is the formation constant. 


Using the formation constants of the iminodiacetic 


given by k_5/K 


acid (IDA) complexes (71) as approximate values for 
tridentate coordination to monoprotonated 1,3-PDTA 
and assuming the analogous formation constants with 
monoprotonated EDDDA to lie between those for IDA 
and iminodipropionic acid, the predicted value for 


the dissociation rate constant of the cadmium inter- 
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mediate with 1,3-PDTA is 1.8 x 10° sec as while that 


with EDDDA is between 1.8 x 10 4 andede 1x0 6 seers 


The predicted value for the dissociation rate constant 


of the zinc intermediate with 1,3-PDTA is 3.7 x 107 
ye while that with EDDDA is between 3.7 x 10° and 
4.5 x 104 Seto The formation constant of the 


lead complex of iminodipropionic acid is not available. 
The reaction by which the proton leaves the 
nitrogen to permit formation of the second metal- 
nitrogen bond could involve either direct transfer 
to solvent water or migration to a neighbouring 
carboxylate group. The kinetics of these two 
reactions for the nitrogen-protonated isomer of 
neutral glycine have been studied by Sheinblatt and 
Gutowsky (72), who report a pseudo first order rate 
constant .of.5 sec | for the transfer of the 
nitrogen-bonded proton to water and a rate constant 
OE nll 3 accu. foremioraction ot the proton to the 
carboxylate group. The rate constants for proton 
transfer in the nitrogen-protonated intermediates 
of 1,3-PDTA and EDDDA would probably be somewhat 
larger than either of those for glycine because of 
charge-charge interaction between the metal ion 
and the proton and the additional carboxylate group 


to which proton. transfer can occur. 
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The preceding considerations suggest that the 
rates of dissociation of the metal ions from the 
EDDDA intermediates are larger than from the cor- 
responding 1,3-PDTA intermediates, and that they may 
be comparable to the rate of proton transfer from 
the nitrogen atom. If so, there is a greater chance 
that EDDDA intermediates will dissociate before 
proton transfer and formation of the second metal- 
nitrogen bond can take place, and consequently k_»5 
would be less for the EDDDA complexes. The preceding 
considerations also suggest that the rate of reaction 
with cadmium should be decreased more than with zinc, 
as observed experimentally. If these conclusions 
are correct, a further decrease would be predicted 
for the formation reactions of the monoprotonated 
form of ethylenediaminetetrapropionic acid. Pearson 
and de Wit (69) found this to be the case for the 


zinc complex. 


Second Order Ligand Displacement Reactions 

The reactions represented by Equations VI to 
VIII are second order ligand displacement reactions. 
The mechanism (7,24,73) for similar EDTA type ligand 
exchange reactions has been demonstrated to involve 


the bound ligand being successively displaced by 
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the incoming ligand with intermediates present in 
which both ligands are simultaneously bound in 
varying degrees to the metal ion. The mechanism 
for the ligand displacement reactions studied in 
this chapter is probably similar. Comparison of 


the rate constants k beewand koe ein TaplelVii for 


SAL114 b 
the cadmium and lead complexes of 1,3-PDTA shows 
a decrease in the values with increasing degree 
of protonation of the incoming ligand, indicating 
that proton transfer might be a rate-determining 


step in these reactions. Similar results have 


been found for other systems (7,39,40,44,73). 
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CHAPTER IV 


EXPERIMENTAL 


A. CHEMICALS 


1,3-PDTA was synthesized by the method of Weyh 
and Hamm (74). The following analytical data were 
obtained for the ligand: Calculated for Cy Hy gN 509: 
es eas ON yee) SROUNU, (>) 22.97 Hy Os Le 
N, 9.1. The 1,4-BDTA was prepared by Dr. George 
Blakney by a procedure analogous to that used for 
1,3-PDTA. EDDDA (Lamont Laboratories) was used 


without further purification. No impurities could 


be detected in the pmr spectra of solutions of each 


of these ligands. Reagent grade metal nitrate 
Salts were used as received. Tetramethylammonium 
(TMA) hydroxide (Eastman Organic Chemicals), tertiary- 


butanol (t-butanol) and 1,4-dioxane (Fisher Scientific) 


were used as received. 


B. SOLUTIONS 


A stock solution of TMA nitrate, used as a 
reference for some of the chemical shift measure- 
ments, was prepared by titration of a 25% aqueous 
solution of TMA hydroxide with nitric acid to a 


neutral pH. 
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The solutions used in the pmr measurements were 
prepared in either triply distilled water or in 
D,0 from the appropriate amounts of metal salt and 
ligand at concentrations varying between 0.10 and 
0.20 M. For making chemical shift measurements 
TMA or t-butanol was added at concentrations of 
0.01 M and 0.03 M respectively to serve as a reference. 
To avoid dilution, the pH or pD was adjusted with 


concentrated HNO KOH or KOD and samples of about 


3! 
0.5 ml were withdrawn at appropriate pH or pD values. 
When kinetic measurements were being made at constant 
pH as a function of free ligand (or metal nitrate) 
concentration, the requisite amount of ligand (or 
metal nitrate) was added, the pH adjusted and a 
sample withdrawn. More ligand was added to the 
solution and the same procedure followed for all 
samples at this pH. The concentrations were 
corrected for the decreasing solution volume. Since 
high concentrations are necessary in the pmr experi- 


ment, no attempt was made to control the ionic 


strength. 


C. pH MEASUREMENTS 


All pH measurements were made at 25°C with an 
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Orion Model 801 pH meter equipped with either a 
standard glass electrode and a fibre-junction 
saturated calomel reference electrode or a micro- 
combination electrode. Saturated potassium acid 
tartrate and 0.01 M sodium tetraborate solutions, 
pH 3.56 and 9.18 at 25°, were used to standardize 
the meter. For D0 solutions, the meter readings 


were converted to pD values using the expression 


of Glascoe and Long (75) 


pD = pH meter reading + 0.40. 


D. PMR MEASUREMENTS 


The proton spectra were obtained at 60 MHz on 
a Varian A-60D, at 100 MHz on a Varian HA-100 
and at 220 MHz on a Varian HR-220 spectrometer. 
These instruments were equipped with Varian variable 
temperature controllers. The temperature of the 
probe was determined using the standard methanol 
and ethylene glycol solutions and the plot of 
chemical shift vs. temperature prepared by Varian. 
At 60 MHz, spectra were recorded at sweep rates 
of 0.1 Hz/sec for both the chemical shift and line- 
shape measurements. The chemical shifts are reported 


in ppm relative to the central resonance of the TMA 
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triplet which is 3.17 ppm downfield from the methyl 
resonance of sodium 3-(trimethylsilyl)-l-propane 
sulfonic acid. For linewidth measurements, each 
spectrum was recorded at least four times and the 
linewidths at half height of the resonances in the 
individual spectra were averaged. 

The HA-100 spectrometer was operated in the 
frequency sweep mode with the methyl resonance of 
internal t-butanol providing the lock signal. Inter- 
ference from the water resonance necessitated 
the use of D,O as the solvent. Each spectrum was 
recorded at least four times at sweep rates of 
0.1 Hz/sec for the linewidth measurements. Fre- 
quencies were measured from the differences between 
the manual and sweep oscillator frequencies. 

At 220 MHz the chemical shift differences and 


linewidths were measured directly from the spectra. 


E. KINETIC APPLICATIONS OF NUCLEAR MAGNETIC RESONANCE 


The basic theories used in the study of chemical 
exchange by nuclear magnetic resonance (nmr) have 
been eee extensively in several monographs and 
review articles (50,76-79). Two of these theories, 


the modified Bloch approach and density matrix theory, 
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have been utilized in the present work to extract 
kinetic data from the pmr spectra. Only a brief 
qualitative description of these two approaches 
will be given. 

Bloch (80) obtained a set of equations to describe 
the variation of the components of the total nuclear 
Magnetic moment per unit volume for a collection 
of nuclei with non-zero spins in the presence of a 
varying magnetic field. The simplest chemical exchange 
system which may be studied by nmr is one for which the 
nuclei may interchange between two different chemical 
environments, A and B with different Larmour fre- 
quencies. For this case, there will be two independent 
aetna ke moments for environments A and B which can 
be described by the Bloch equations. Addition of terms 
to these equations to account for the transfer of 
magnetization between the two sites leads to the 
modified Bloch equations (81). From these equations 
a general expression may be derived describing the 
intensity of the nmr absorption at a given frequency 
when there is exchange between the sites (82,83). 

This general equation may be simplified under certain 
limiting conditions. 

One such condition is that of slow exchange for 


which the spectrum consists of two separate resonances 
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at frequencies 5 andes "atone a £Or chilis situa- 


B 


tion the lifetimes TH and Tp of the nuclei in each 


of the sites before exchange to the other site are 


large compared to the inverse of the separation of 


-] 
Ray oh (5, 53) . These 


lifetimes may be determined from the extent of broad- 


the resonances; that is, 


ening of each of the resonances by an equation of the 


EON fo); 
i. 2 
apeea ts naw ek aaa 
' 
where Wiy2 A is the linewidth at half height of the 
; 44 


exchange — broadened resonance of nuclei at site A 

and Wive A is the linewidth of the same resonance in 
’ 

the absence of exchange. An analogous equation 


may be written for the lifetime t The natural 


5° 
linewidth W is related to the effective spin- 
1/2hA 


spin relaxation time T by the relation 1/T = 


2,A 


7™W - The effective spin-spin relaxation time 


2jA 


takes intg account the actual spin-spin relaxation 
‘time and broadening due to magnetic field inhomo- 
geneities. 

For the condition of intermediate exchange, the 
lifetimes are of the order of the inverse of the 
separation of the resonances and the spectrum consists 


of the transition from two resonances to one. When the 
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lifetimes are small compared to the inverse of the 
separation of the resonances, rapid exchange occurs 
and the spectrum consists of one resonance at an 
intermediate frequency between by and Sp: 

Most of the ligand exchange systems studied in 
Chapter III are ones for which Equation XVIII is 
applicable. The two sites are free ligand and com- 
plexed ligand. 

The modified Bloch equation approach may be 
extended to the many-site exchange problem by the 
addition of appropriate exchange terms and the 
derivation of an expression describing the intensity 
of the nmr absorption as a function of the frequency. 
The lifetimes in the individual sites may then be 
obtained by computing theoretical spectra and ave ee 
them to the experimental spectra. Such an approach 
was used for the Pb-1,3-PDTA exchange system which 
was analyzed as a four-site exchange problem. 

The exchange problem becomes much more complex 
when there is spin-spin coupling between the nuclei 
at the different sites. For such systems the 
modified Bloch approach is no longer valid; however, 
a quantum mechanical method based on the use of 
density matrices is applicable. In this theoretical 


treatment (78,79), exchange and relaxation matrix 
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terms are added to the basic density matrix equation 
using as the basis set the spin product functions 
for the spin system in question. A final lineshape 
expression for the nmr absorption spectrum can then 
be derived. For an AB spin system, the lifetimes 

tT, Or Tt, may be obtained by computing theoretical 
spectra using the above lineshape expression and 
matching them to the experimental spectra. Such an 


approach was used in Chapter II for determining the 


lifetimes before AB interchange. 
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THE BINDING OF METAL IONS BY GLUTATHIONE 
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CHAPTER V 


INTRODUCTION 


The binding of metal ions by molecules such as 
peptides and proteins is of fundamental interest in 
view of the importance of metal ions in biological 
systems (1-4). For example, the catalytic function 
of many enzymes is dependent upon the presence of metal 
ions which may be coordinated to certain functional 
groups of the enzyme or substrate. In an enzyme 
many groups are available as potential coordination 
Sites, but only certain ones are actually metal- 
coordinated in an enzymatic reaction. 

Metal-peptide complexes are potential models for 
metal-enzyme interactions since peptides can be syn- 
thesized which contain the functional groups involved 
at the active sites of enzymes. Moreover, peptides 
are relatively small molecules making them more 
amenable to study by a variety of techniques. Freeman 
has summarized X-ray crystallographic data for many 
metal-peptide complexes (5). Even though these data 
refer to solid state structures they may be useful for 
at least suggesting possible metal binding sites and 
conformations of the complexes in solution. A great 


deal of research has also been devoted to the solution 
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chemistry of metal-peptide complexes. Several of these 
studies will be reviewed briefly below. 

Representative of the early work in this area are 
the studies in which the formation constants of the 
cadmium, zinc, and lead complexes of a series of | 
polyglycine peptides were determined (6-8) using 
potentiometric measurements. By comparing the formation 
constant data, it was proposed that the binding sites 
in these cOmpiexes consist of the terminal amino group 
and the adjacent peptide linkage, but it was not possible 
to establish any carboxyl coordination. The copper 
and nickel complexes of the polyglycine peptides have 
also been extensively investigated (9-12,18) using 
pH titration, infrared and spectrophotometric methods 
in order to determine formation constants. Ionization 
of the peptide protons and subsequent metal binding 
to the negative nitrogen atoms was proposed by some 
of the workers for these complexes. pH titration 
and spectrophotometric studies have also been carried 
out on the copper complex of s-alanyl-L-histidine 
(carnosine) (9,13,14). From this work it was suggested 
that ionization of the peptide proton takes place 
with subsequent copper binding to the negatively 
charged peptide nitrogen, but there was lack of agree- 


ment as to the binding sites and the structure of 
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the complex. 

The techniques used in the above examples have 
the disadvantage that they measure macroscopic proper- 
ties and hence, do not provide definitive information 
at the molecular level. Pmr spectroscopy has proven 
useful for elucidating the binding of metal ions by 
simple peptides at the molecular level. For example 
the functional groups of polyglycine peptides involved 
in coordination to diamagnetic cadmium, zinc, lead 
and nickel (15-17) have been established from changes 
in the chemical shifts of carbon-bonded protons 
close to the binding site. The groups of polyglycine 
peptides involved in binding to paramagnetic copper 
and nickel (16,17) and of carnosine with copper (38) 
have been identified from the dependence of the line- 
widths of the pmr resonances on their proximity to 
the binding site. This application of pmr is limited, 
however, to those peptides whose spectra are relatively 
simple because of the need for distinct, well resolved 
resonances for monitoring interactions at the potential 
binding sites. Spectra which consist of overlapping 
peaks and complex multiplet patterns are of limited 
use in metal binding studies. 

Carbon-13 magnetic resonance (cmr) spectroscopy 


is potentially more useful than pmr for the study of 
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metal-peptide interactions. First of all, proton- 
decoupled cmr spectra are comprised of single 
resonances for each of the nonequivalent carbon atoms 
of peptides. Secondly, the range of chemical shifts 
in cmr is at least an order of magnitude greater 
than in pmr, so that often a resonance can be re- 
solved for each of the carbon atoms. In addition, 
carbon atoms are bonded directly to the possible 
binding sites in peptides making cmr potentially 
more sensitive as a probe for studying binding at 
the molecular level. 

li cnepter Vil, the Lesutrts ol “a cmr”“investi- 
Cael Or The Dilding OL cadmium, Zinc, Lead and 
mercury by the tripeptide y-L-glutamyl-L-cysteinyl- 
. glycine: (glutathione) are presented and discussed 


(23). This peptide was chosen to evaluate the 
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potential of cmr for the elucidation of metal binding 
by peptides because there is a lack of agreement as 
to which of the six possible coordination sites are 


involved in binding to these metal ions (19-22). 
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CHAPTER VI 
THE BINDING OF ZINC, CADMIUM, LEAD AND MERCURY 


BY GLUTATHIONE 


The tripeptide glutathione, which is widely 
distributed in nature, contains the following six 
potential sites for metal ion coordination: the 
sulfhydryl group, the amino group, the two carboxyl 
groups and the two peptide linkages. Formation 
constants for the zinc, cadmium and lead complexes 
of glutathione have been determined previously 
Mean. Cheap. titrations method. (19-22), but_the 
results are different in each case because of the 
lack of agreement as to which of the sites are 
involved in metal ion coordination. The validity of 
microscopic formation constants determined from pH 
titration data is questionable when the binding sites 
are not known because the method is based on the com- 
petition between the metal ion and the proton for 
the site. Furthermore,amicroscopic Lonization 
constants must be used in the calculations and 
only ingone instance (21) was this done. 

Li and coworkers (19,20) obtained formation 
constants for the zinc, cadmium and lead complexes 
of glutathione by assuming that all ionizable 


protons are removed from the peptide when it binds 
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to these metal ions. They suggested that in these 
complexes the metal ion is bound to the ionized 
sulfhydryl and amino groups. Martin and Edsall (21) 
on the other hand, considered that divalent metal 
ions could be bound to glutathione either through 
the amino and carboxyl groups of the glutamyl residue 
or to the sulfhydryl group of the cysteinyl residue 
with possible chelate formation through a peptide 
bond. By assuming that sulfur binding could be 
excluded in S-methylglutathione and that bonding to 
the glutamyl residue would be equivalent in gluta- 
thione and its S-methyl derivative, they assigned 
stability constant values to the sulfhydryl bound 
and to the glutamyl bound complexes of zinc. At 

pH < 8 these workers postulated that zinc binding 
could take place at a peptide bond, but that ioniza- 
tion of a peptide proton was unlikely. Perrin and 
Watt (22) considered a number of different zinc and 


cadmium complexes of glutathione of the type ML , 


MLo MHL, ML(HL)>, and calculated formation con- 
stants for each of these species. These workers 


proposed structures involving simultaneous coordina- 
tion of sulfhydryl and glutamyl groups to the same 
metal ion, and suggested that ionization of the 


peptide protons can occur on complex formation. 
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Because of the differing results of the previous 
studies of the metal complexes of glutathione, this 
particular peptide was chosen to evaluate the potential 
of cmr for the elucidation of metal binding by pep- 
tides. The results of a cmr investigation of the 
binding of zinc, cadmium, lead and mercury by gluta- 
thione are reported in this chapter. Binding to 
mercury has been studied previously by the polaro- 


graphic method (24), and binding to CH Hg” has been 


3 
studied by pmr (25). In addition, the acid-base 
chemistry of the four acidic groups of glutathione, 
which is necessary for quantitative studies of metal 


binding, has already been characterized at the 


molecular level (25,26). 


A. RESULTS 
Glutathione 

Cmr spectra of glutathione were obtained over 
a pD range WAtoes 3 sror ab. 30M solution om the, trig 
peptide in D.0 using 1,4-dioxane as an internal 
reference. The spectrum obtained at pD 7.37 is the 
one designated (a) in Figure 15. The carbon atoms 
are identified by Glu, Cys or Gly to specify the 
amino acid residue glutamyl, cysteinyl or glycyl in 


which the carbon is located and Cur & Cr CONH or 
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Figure 15: Sd magnetic resonance spectra of gluta- 


thione. 


(a) D0 solution containing 0.30 M gluta- 


thiene-at pd-7.374 


(b) D0 solution containing 0.30 M gluta- 


thione and 0.15 M Zn (NO,) . at pD 7.84. 
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COOH to indicate the particular carbon of that 
residue. The assignments of the carbon resonances 
were made by comparison with the chemical shifts of 
the constituent amino acids (27,28), and from the 
chemical shift dependence of each of the resonances 
on the pD. These assignments are in agreement with 
those reported by Jung and coworkers who have 
previously given the cmr spectrum of glutathione 
(29), and have presented chemical shift vs. pH 


curves for each of the carbon atoms (30). 


Zinc-glutathione 

Addition of zinc to a solution of glutathione 
causes the chemical shifts of selected resonances 
to change, as illustrated by spectrum (b) en Figure 15. 
This spectrum was obtained from a solution 0.30 M 
in glutathione and 0.15 M in Zn (NO,) eth e! D,O at 
DORs Bawa enol CaLes. that. the ac chemical shifts are 
sensitive to metal binding. The resonances were 
assigned by following their pD dependence from low 
pD values where little metal binding takes place. 
The chemical shifts of the three cysteinyl carbon 
resonances are shown as a function of pD for a zinc 
to glutathione ratio of T:2 by the solid curves ian 


Figure 16, while the chemical shifts of selected 


carbon atoms of the glutamyl and glycyl residues 
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CHEMICAL SHIFT (ppm from dioxane} 
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Figure 16: 


ied 








pD dependence of the chemical shifts of the 
cysteinyl carbons for a D.0 solution con- 
taining 0.30 M glutathione and 0.15 M 

Zn (NO) 5- The dashed curves represent the 
chemical whee behavior of the same carbon 


atoms for a solution of glutathione con- 


taining no complexing metal ion. 
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are shown for the same solution conditions in Figure 
17. The discontinuity in these solid curves indicates 
the pD region over which precipitation occurred. 
The dashed curves in Figures 16 and 17 represent the 
chemical shift vs. pD behaviour of the indicated 
carbon resonances in the absence of coordinating 
metal ion. 

The lack of any large dependence of the chemical 
shifts plotted in Figures 16 and 17 on the presence 
Of zanceat pD < 32 indicates no strong coordination 
in this pD range. A small amount of binding to both 
carboxylic acid groups, which previous studies (25) 
have shown to be partially ionized at pD > 0, is 
indicated by the small chemical shift differences 
observed for the Gly-C., Gly-COOH, Glu-C and Glu-COOH 
carbons. A chemical shift difference for a given 
carbon atom refers to the difference between the solid 
and dotted. curves in Figures 16 and 17 and in sub- 
sequent figures, and is assumed to be indicative 
of some degree of metal binding to a specific site 
near the given carbon atom. Some carboxylic acid 
group coordination is consistent with the known 
binding of zinc by the carboxylic acid groups of 
acetylglycine and the C-terminal end of polyglycine 


peptides (15,31). 
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CHEMICAL SHIFT (ppm from dioxane) 


Figure 17: 





Glu - CONH 


pD dependence of the chemical shifts of 
selected glutamyl and glycyl carbons 
pal a hears) D0 solution containing 0.30 M 


glutathione and 0.15 M Zn (NO3).,. 
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Between pD 3.0 and 13.2, the large chemical 
shift differences exhibited by all three cysteinyl 
carbon resonances indicate that zinc is bound to some 
extent to the sulfhydryl group and possibly to the 
peptide linkage between the cysteinyl and glycyl 
residues. The lack of any chemical shift difference 
for the Glu-CONH carbon at pD < 10.5 indicates no 
binding to the peptide linkage between the glutamyl 
and cysteinyl residues. The Glu-COOH resonance 
indicates a small amount of binding to the glutamyl 
residue up to pD 6, presumably involving only the 
carboxyl group. An increase in the binding to the 
glutamyl residue at pD > 6, indicated by the chemical 
shift differences for the Glu-COOH, =C RY -C., and ace 
carbon resonances, probably involves simultaneous 
Zinc coordination to the amino and carboxyl groups 
by analogy with the coordination of zinc by glycine. 
The small chemical shift differences observed for 
the Gly-C. and Gly-COOH resonances over the pD range 
3 to 9 are consistent with a small amount of binding 
to the glycyl carboxylic acid group. 

The rather large downfield shifts exhibited 
by the Cys-CONH, Cys-C., 
resonances at pD > 10.5 indicate a change in the 


Glu-CONH and Steere carbon 


nature of the binding to the peptide linkages. If 
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binding does take place at the peptide linkages at 
pD < 10.5, it is presumably to the carbonyl oxygen 
of the neutral peptide linkage, since it is known 
to be the more basic site (5,32-34). The downfield 
shift behaviour suggests that, in the presence of 
zinc, ionization of one or both of the peptide 
protons might occur at pD > 10.5 with subsequent 
binding of zinc to a negatively charged nitrogen 
atom. Separate ee ese for the free amino acids 
of which glutathione is comprised were not observed, 
indicating that hydrolysis of the peptide bonds had 
not occurred under the given solution conditions. 
The specificity of the zinc binding in the JpD 
range 3 to 6 was investigated further by monitoring 
the chemical shifts as a function of the zinc to 
glutathione ratio for an H,0 solution at a constant 
pH of 5.51; the results of this study-are "presented 
in Figure 18. It can be seen that the chemical 
SHLLC of che Cys-C, carbon varies continuously as 
the ratio is increased from 0.0 to 0.5 and then 
levels out at a constant value of 38.2 ppm; the total 
chemical shift variation is 2.7 ppm. This indicates 
that at a ratio of 0.5 and a pH 62 .5.517) (fhe main 
species in solution is a complex consisting of one 


zinc and two glutathione molecules. The chemical shift 
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Figure 18: The chemical shifts of selected carbons 
as a function of the zinc to glutathione 


ratio for an aqueous solution at pH 5.51. 
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change of only 0.3 ppm for the Glu-COOH carbon up 

to a ratio of 0.5 suggests that only a small amount 
of binding takes place at the glutamyl end. For 
comparison, ionization of the glutamyl ammonium 
proton results in a change of 8.4 ppm for this carbon 
resonance. On-Changing the ratio from 0%.5--to-0.7, 
however, when the sulfhydryl group is presumably com- 
pletely complexed, the additional chemical shift 
change of 0.6 ppm for the Glu-COOH carbon indicates 
that the glutamyl end is becoming more important as 

@ binding site. The™chemical shift variation of. only 
1 ppm for. the Gly-C carbon in Figure 18 indicates 
very weak glycyl carboxyl binding at the given pH. 


Precipitation occurred at ratios greater than 0.7. 


Cadmium-glutathione 

The chemical shifts of the three cysteinyl 
carbons for a cadmium to glutathione ratio of 1:2 
are shown as a function of pD in Figure 19 by the 
solid curves. The corresponding curves for selected 
carbon atoms of the glutamyl and glycyl residues are 
shown in Figure 20. The dotted curves in these 
A cars represent the data points for a glutathione 
solution containing no complexing metal ion. At 
pD < 2, evidence for a small amount of cadmium 


binding to the two carboxyl groups is given by the 
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CHEMICAL SHIFT (ppm from dioxane) 





Figure 19: pD dependence of the chemical shifts of 


the cysteinyl carbons for a D.O solution 


2 
containing 0.30 M glutathione and 0.15 M 


Cd (NO3).- 
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20: pD dependence of the chemical shifts 
of selected glutamyl and glycyl carbons 


for a D,O solution containing 0.30 M 
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small chemical shift differences observed for the 
Gly-C., Gly-COOH, Glu-C, and Glu-COOH carbons. As 
in the zinc system, this type of binding is consistent 
with the known binding of cadmium by the carboxylic 
acid groups in certain peptides (31). 

Between pD 2.0 and heer large chemical 
shift differences observed for the three cysteinyl 
carbon resonances indicate that cadmium is bound 
to the sulfhydryl group and possibly to the peptide 
linkage between the cysteinyl and glycyl residues. 
No detectable binding to the other peptide linkage 
is indicated by the fact that the chemical shift of 
the Glu-CONH carbon resonance is not changed by the 
presence of cadmium. The chemical shift curves 
for the Glu-COOH carbon suggest that a small amount 
of glutamyl binding takes place up to pD 7, probably 
involving the carboxyl group. An increase in glutamyl 
binding at pD > 7, indicated by the chemical shift 


B 


resonances, most likely involves simultaneous 


behaviour of the Glu-COOH, “Cor -—C and me) carbon 


coordination of cadmium to both the glutamyl amino 

and carboxyl groups. Finally, the small chemical 
shift differences observed for the Gly-C, and Gly-COOH 
carbons over the pD range 2 to 9 indicates weak 


glycyl carboxyl coordination. 
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The above results indicate that, in the pD range 
2 to 7, cadmium binding occurs almost exclusively at 
the cysteinyl residue of glutathione. This specifi- 
city was further investigated by monitoring the chemical 
shifts as a function of the cadmium to glutathione 
ratio for an H,0 SOlLULLOn ote wonetant oh of '6..59- 
The chemical shifts of selected carbon resonances 
are plotted as a function of the ratio in Figure 21. 
The cysteinyl residue is the principal coordination 
site up to a ratio of 0.5, as evidenced by the: chemical 
Shift variations of the cysteinyl carbons. In par- 


tichley~s tue chemical shirt of the Cys-C. carbon «changes 


B 
by 3.2 ppm and then levels out above a ratio of 0.5. 
The chemical shift of the Glu-COOH carbon changes by 
only s0.> ppm Dp co thtesratio,indicating very. ttle 
glutamyl binding. However, an increase in this inter- 
action takes place at ratios greater than 0.5, when 
presumably the sulfhydryl group is fully coordinated. 
Finally, the chemical shift change of only 0.2 ppm 

for the Gly-C resonance indicates weak binding to 


the glycyl carboxyl group under these conditions. 


Precipitation occurred at ratios greater than 0.7. 
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CHEMICAL SHIFT (ppm from dioxane) 


Glu-COOH 





0.0 0.2 0.4 0.6 
[(Cd2*]/ [GLUTATHIONE] 


Figure 21: The chemical shifts of selected. carbons as 
a function of the cadmium to glutathione 


ratio. for, ansagueous, Solution. at pH 6559. 
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Lead-glutathione 


The chemical shift curves of the three cysteinyl, 
the Gly-C and the Glu-COOH carbons are shown as a 
function of pD in Figure 22 for a lead to glutathione 
Katiovol Ure s precipitation from “pp 2t oto 5.4 and 
above pD 12 limited the range over which binding could 
be studied. 

At pD < 2.3, weak coordination to the two carboxyl 
groups is indicated by the small chemical shift dif- 
ferences for the Gly-C and Glu-COOH carbons. Between 
pD 5.4 and 12, the chemical shift curves provide 
evidence for lead binding to only the cysteinyl and 
glycyl residues. The large differences between the 
solid and dotted curves for the three cysteinyl carbon 
resonances indicate lead coordination to the sulfhydryl 
group and possibly the peptide linkage between the 
cysteinyl and glycyl residues. Evidence for weak 
glycyl carboxyl binding between pD 5.4 and 9.0 is also 
provided by the chemical shift curves for the Gly-c, 
and Gly-COOH carbons. The extent of this binding 
decreases at pD > 9 presumably due to carboxyl displace- 
ment by hydroxide ions with the formation of lead- 
hydroxy-glutathione mixed complexes. The lack of any 
change in the chemical shifts due to the presence of 


lead for all the glutamyl carbon resonances is evidence 
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Figure 22: pD dependence of the chemical shifts 


of selected carbons for a D0 solution 


containing 0.30 M glutathione and 


0.15 M Pb(NO,),. 
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that the glutamyl functional groups do not act as 


binding sites in the pD range 5.4 to 12. 


Mercury-glutathione 

The chemical shift curves for the three cysteinyl 
carbons as a function of pD, shown in Figure 23 for a 
mercury to glutathione ratio of 1:2, are seen to be 
nearly constant over the entire accessible pD range. 
The chemical shifts for all the other carbon atoms of 
glutathione are identical with those obtained from 
solutions containing no mercury ions. Therefore, the 
binding of mercury appears to take place exclusively 
to the sulfhydryl group at a mercury to glutathione 
ratio of 0.5. The results of further experiments, in 
which the ratio was \wariedtirom 0.0 to 0.5 for an H,O 
solution at a constant pH of 5.01 are presented in 
Figure 24. These plots confirm the specificity of the 
sulfhydryl group for mercury. Precipitation prevented 
investigations at mercury to glutathione ratios 


greater than 0.5. 


B. DISCUSSION 


The cmr results show that all of the metal ions 
studied in this work bind to the sulfhydryl group of 


glutathione at metal to peptide ratios of 1:2 and at 
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Figure 23: pD dependence of the chemical shifts of 


the cysteinyl carbons for a D0 solution 


containing 0.30 M glutathione and 0.15 M 


Hg (NO3) 5 





ee ee ee ae Oe ee ee Se an 
; - _ 
al peo ~ 


dq 7 SS we Ree the - 
lo adline Incimed> st to sonebasqeh Gq sta a1 : 
notsuloe G.4 s 103 etodis5 Iyntedeyo sat 


seh AREAS) how anciittegate 1°0C 0 acaba 
ee ha eh! ; 7 






é 5 
— 
_ 





oy we 
- 


137 





> 40 

| = 

oO 

x 

oO 

GS 38 

6 < Cys-Ce 

E 

oO: 

igs 2 

= CATS 6 
LU. 19 eee 
ae 

2) 

q -106 Cys-CONH 
Oo 

Ss 

oui Nhe Glu-COOH 
©O 


—NyO Gly- COOH 


0.0 0.1 0.2 0.3 0.4 0.5 
[Hg*"]/ [GLUTA THIONE] 


Figure 24: The chemical shifts of selected carbons 
as a function of the mercury to gluta- 
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pD > 3. Also, all of the metal ions except mercury 
coordinate rather weakly to the glycyl and glutamyl 
carboxyl groups, while only zinc and cadmium bind to 
the glutamyl amino group. The extent of binding to 
the various sites is dependent upon the solution pD, 
except in the case of mercury which is bound strongly 
to the sulfhydryl group over the entire accessible pD 


range. 


Zinc and Cadmium Binding to the Sulfhydryl and Glutamyl 
Groups 

Because of the similarity in the binding of cadmium 
and zinc, the results for these systems will be dis- 
cussed together. The strongest binding of these two 
metal ions occurs to the sulfhydryl group, as evidenced 
by the fact that the metal ions bind to this group at 
a much lower pD than to the glutamyl amino group, 
even though the PK, of the sulfhydryl group is only 
Oe 2eiiiits less than that. of thie amino. group 6(25). » Th} 
chemical shifts of the Cys-C,, carbon indicate that, 
foros metal to oluitathione, batio ote 0.5, al) woeistir 
hydryl groups are complexed at pD 6 in the zinc system 
and at pD 7 in the cadmium system. The results of 


experiments in which the metal to glutathione ratio was 


varied at pH 5.51 for the zinc system and at pH 6.59 
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for the cadmium system also indicate that binding 
occurs almost exclusively to the sulfhydryl group at 
these. pH, values..up to: a ratio of..0.5... Thus, a 1:2 
complex such as that shown below is the major species 


in solution for these conditions. 


—o © 
2 CCHCH.CH.,CONHCHCONHCH -CO 


& 
NH3 sie 
5 
| 
M 
| 
5 
| 
CH, 
© Conch. GH paedas meme CO © 
2 Moe, ee eo 
© 
NH, 


This 1:2 complex is analogous to the complex 
M(HL) 5, where vee is eee or ie and cnr represents 
monoprotonated glutathione, proposed by Perrin and 
Watt (22). These workers were unable to assign the 
binding sites due to the uncertainty in the PK, 
values of the amino and sulfhydryl groups. The amino 
group is still protonated as evidenced by the chemical 
SHitt curves ror the Glu-COOHu, ae and “Ce carbons. Li 
and coworkers (19,20) interpreted pH titration data 
for the binding of cadmium and zinc by glutathione in 


terms of complexes in which the two carboxyl, the amino 
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and the sulfhydryl groups are all ionized. The present 
chemical shift data show this not to be the case, but 
rather that several different complexes may form and in 
some, those functional groups which are not coordinated 
are protonated. 

It is possible that, in the zinc and cadmium 
complexes discussed above, binding might also occur 
to the peptide linkage between the cysteinyl and glycyl 
residues. Of the two sites in this linkage, protonated 
nitrogen and carbonyl oxygen, the oxygen is the more 
basic (5,32-34) and hence a more probable binding site. 
Such binding occurring simultaneously with sulfhydryl 
group coordination would result in a six-membered 
chelate ring. The cysteinyl carbon chemical shift 
data for the zinc system at pD 7 are consistent with 
carbonyl oxygen binding. 

As shown by the chemical shift vs. pD curves, 
observable metal binding to the glutamyl amino group 
takes place at pD > 6 for the zinc system and at pD > 7 
for the cadmium system. The extent of this binding is 
pD dependent, as illustrated by the results in Table 
XVII which were calculated from the chemical shift data 
for the Glu-COOH and Glu-C, carbon resonances using 


B 


the equation, 
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Table XVII 


Fraction of Glutamyl Groups of Glutathione 


Complexed by ane and ca2t as* a’ Function, of? pD 3 
gn2t ca2t 

pD Glu-cooH® Glu-C, SS lecin-coon Glu-C, 
Ge (a Pa 0 ie 0.00 0; 03 0.00 
rir 0.10 Vee 0.05 0.00 
Ze 0.34 0.29 0.14 0.06 
Big 0.49 0.41 0.21 0.14 
9. 0.58 Oras eos 0236 

10. Us oe tps) Ooy Ue 

ee Use 075.6 epee ieee 

(a) Data obtained from solutions of metal to glutathione 

YAcCio, OL 132. 
(b) Glu-COOH carbon chemical shift data. 


(c) 


SEE carbon chemical shift data. 
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where f is the fraction of glutamyl groups complexed, 
Sons is the observed chemical shift of a glutamyl carbon 
in the presence of complexing metal ion, b- is the 
chemical shift in the absence of metal ion and Pe is the 
chemical shift of glutamyl complexed glutathione. The 
value of Se. used in these calculations was obtained from 
the point where the solid and dashed curves cross in 
Pisures sjeand 20.(31). .ithe Sons values were obtained 
from the chemical shift vs. pD curves for solutions 

of a metal to glutathione ratio of 1:2, and those for 

5 - from the corresponding curves for a solution of 
glutathione. The fractions shown in Table XVII reveal 
that cadmium is bound less strongly to the glutamyl end 
CragtuloatLbionesrhan i ss Zinc. 

A change in the nature of the zinc binding to the 
evetoinyvlsaresiduelsin the pD range 6, t0;.10. 54 1s,andicated 
by the chemical shift data for the cysteinyl carbon 
resonances shown in Figure 16. Simultaneous coordination 
to the sulfhydryl group and the peptide carbonyl oxygen, 
as suggested above, would be expected to have a different 
effect on all three cysteinyl carbons than coordination 
to the sulfhydryl group alone. Therefore, it is sug- 
gested that at pD < 6, simultaneous sulfhydryl-peptide 
oxygen binding takes place, and at pD > 6 the amount of 


binding to the peptide oxygen decreases. Presumably 
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the decreased peptide binding is due to glutamyl group 
competition for the coordination sites on the zinc. 
The smaller affinity of cadmium for the glutamyl group 
probably causes a smaller decrease in peptide binding 
and hence, the effect on the cysteinyl carbon chemical 
shifts is less than in the zinc system. 

At pD > 6 for the zinc-glutathione system and pD > 7 
for the cadmium system, the principal binding sites are 
the glutamyl residue and the cysteinyl residue. Glutamyl 
binding most likely involves simultaneous metal co- 
ordination to the amino and carboxyl groups by analogy 
with the coordination of these two metal ions by glycine. 
As indicated above, cysteinyl binding could involve 
Simultaneous coordination to the sulfhydryl group and 
the oxygen of the peptide linkage between the cysteinyl 
and glycyl residues. In any case, it is unlikely that 
both the cysteinyl and glutamyl residues of a glutathione 
molecule are simultaneously coordinated to the same 
metal ion; this would involve an unstable 10-membered 
ring (35), since the present chemical shift data have 
shown that the peptide linkage between the glutamyl and 
cysteinyl residues is not bonded to cadmium nor to 
zinc at pD < 10.5. The structure proposed by Perrin and 


- 4 + 
Watt (22) for the ML complexes, where ven! is ca? or 


anes and a represents fully ionized glutathione, 


aa 
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involves simultaneous coordination to the glutamyl 
amino and carboxyl groups, the glutamyl peptide linkage 
and the sulfhydryl group. This clearly is not in 


agreement with the L3¢ chemical shift data. In the 


structure proposed for the MLo complexes by these 
workers, the second glutathione molecule was assumed 

to be bonded to the metal only through the glutamyl 
carboxyl and amino groups. More likely structures 

for these complexes are ones in which the cysteinyl 

or glutamyl group of one glutathione molecule and 

the cysteinyl or glutamyl group of another are bound 

to the same metal ion, as considered by Martin and 
Edsall. (21) \inethe,janalysis+of pH titrationidata. 
Another likely structure might be one in which the 
cysteinyl and glutamyl groups of a glutathione molecule 
are bound to different metal ions, resulting in 
polynuclear complexes which have been proposed 
previously for polyglycine peptide complexes (15,16). 
Molecular models indicate that polynuclear complexes in 
which a metal is bonded to the sulfhydryl groups of 

two glutathione ligands and the glutamyl amino and 
carboxyl dentates of two other sulfhydryl-complexed 


glutathione molecules are possible. 
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Ionization of the Peptide Protons 

It is known that, in certain metal-peptide 
complexes, the metal ion promotes the ionization of 
peptide protons with subsequent binding to the negatively 
charged nitrogen atom (9,10,16,18,33,36,37). The 
large downfield shifts observed for the Glu-CONH, 
aa Sy, Cys-CONH, and Cys-C,, carbon resonances for the 
zinc system at pD > 10.5 suggest that zinc might promote 
ionization of the peptide protons of glutathione. After 
such ionization the most basic site is the nitrogen 
atom (5,34) where zinc binding could take place. 
Simultaneous zinc coordination to the sulfhydryl group 
and glutamyl peptide nitrogen would result ina five- 
membered ring and to the cysteinyl peptide nitrogen 
would result in a six-membered ring. Ionization of a 
peptide proton and binding of the negatively charged 
nitrogen atom to a zinc ion presumably have different 
effects on the chemical shifts of the neighbouring 
carbon atoms. The fact that the chemical shift of 
the Cys-C. carbon, which is near both peptide nitrogens, 
experiences no change at pD > 10.5 might be due to 
the combined effects of ionization and binding can- 
celling each other out. The chemical shift data for 
the cadmium, lead and mercury systems provide no 


evidence for peptide proton ionization. 
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Perrin and Watt (22), from pH titration data, 
suggested that ionization of the peptide protons of 
glutathione can occur upon complex formation with both 
zinc and cadmium and that the metal-peptide bond 
interaction is slightly greater in the zinc system. 
However, the PK, of 9.86 which they interpreted as 
corresponding to ionization of a peptide proton from 
a coordinated peptide group in the zn complex, 
is approximately two orders of magnitude too low to 
be consistent with the present chemical shift data. 
Martin and Edsall (21) found that their pH titration 
data for zinc-glutathione did not necessitate the 
postulation of the ionization of peptide protons; 
however their experiments were carried out at pH < 8. 
These same workers have suggested that zinc ion does 
promote peptide ionization in the zinc-glycylhistidine 


complex (13). 


Mercury and Lead Binding to the Sulfhydryl Group 


The binding of mercury and lead to glutathione 
is somewhat different from that of zinc and cadmium; 
neither’-mercury nor lead binds to the glutamyl end 
abupD oc cLorxmetiah ito wpeptide Inatiios wf, QeoaisMercury 
binds exclusively to the sulfhydryl group at ratios up 


to 0.5 over the pD range 0.5 to 12, resulting in 1:2 
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complexes with the protonation state of the carboxyl and 
amino groups dependent on the solution pD. No binding 
to any of the other potential coordination sites is 
indicated by the chemical shift data. This is in agree- 
ment with the results of Stricks and Kolthoff (24) 


who found that in the complex HgL where L represents 


pid 
glutathione in any of its protonation states, mercury 
is firmly bound to glutathione as a mercaptide. 


Since mercury and CH Hg? (255 39) nbandrtostheysulfhydrvyl 


a 
group with such a high degree of specificity and since 
such binding results in large chemical shift changes 


in the Cys-C. and -C, resonances, it may be possible 


B 
to identify these two resonances in the cmr spectra 
of proteins by observing changes in the spectrum as 


: ; : ; + 
the protein is titrated with mercury or CH,Hg at 


3 
constant pH. 

The binding of lead to the sulfhydryl group of 
glutathione is not as strong as that of mercury, as 
evidenced by the decrease in the extent of binding 
atipdisn? andpby7thetfiacththatehosbindinghtp the 
sulfhydryl group is observed at pD < 2.3 for the lead- 
glutathione system. The chemical shift data for the 
Cys-CONH carbon indicate that lead might be simul- 


taneously bound to the sulfhydryl group and to the 


oxygen of the peptide linkage between the cysteinyl and 
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glycyl residues. The chemical shift data also indicate 
that the glutamyl amino group is not a binding site 
for a lead to glutathione ratio of 0.5 over the 


accessible pD range. 


Chemical Shift Differences 

Addition of metal ions to a solution of glutathione 
causes the chemical shifts of selected carbons close 
to the binding site to change. Results of gquantita- 
tive measurements of the magnitudes of the differences 
between -the chemical shifts of selected carbon atoms, 
when a given functional group is complexed and when the 
group is not complexed, are given in Table XVIII. The 
magnitudes of these chemical shift differences are 
denoted by A. 

In the determination of the A values for the 
cysteinyl carbon resonances, the chemical shifts when 
the sulfhydryl group is complexed were obtained from 
the chemical shift vs. pD curves for the metal Eon ee 
solutions at pD 7. The chemical shifts when the 
sulfhydryl group is not complexed were obtained from 
the chemical shift vs. pD curves for glutathione 
solutions at pD 7, where the sulfhydryl group is 
protonated. The order of decreasing 4 values for the 


2+ 2+ 2+ 


+, 
Cys-C, carbon is Hg > PD aad > ans , in agree- 
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Table XVIII 


1 
36 Chemical Shift Changes for Selected Carbon 


Resonances of Glutathione upon Complexation 


A (ppm) 
CARBON gn°* cat Pan Hg** 
Cys-C., ert: Rhee bi 4.4 
Cys-C Hr 8 bes bes 1.0 
Cys-CONH a7 Uns 0.8 Vey 
Glu-COOH 6.0 ae a a 
Glu-C,, 20 (caps a a 
Gly-C. O51 OFZ ae b 


(Ase NOecCOOLCd NAL On. Or Hg?* or Pho? to the glutamyl 
amino group was detected. 
(ao) No. coordination of Hg? to the glycyl carboxyl 


group was detected. 
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ment with the order of decreasing affinity of the 
sulfhydryl grouv of bovine serum albumin for these 
metal ions (40). The data in Table XVIII also 
illustrate that the A values decrease as the number 

of bonds separating the carbon atom from the sulfhydryl 
group increases. The much greater decrease for the 
Cys-C, and Cys-CONH carbons in the mercury system is 
consistent with no binding of this metal ion to the 
oxygen of the peptide linkage between the cysteinyl 

and glycyl residues. 

The A values for the Glu-COOH and Glu-C, carbon 
resonances are presented in Table XVIII for the cadmium 
and zinc systems; no glutamyl amino coordination by 
mercury or lead was detected by the chemical shift 
data. The chemical shifts for the glutamyl complexed 
species were obtained from the points where the solid 
and dashed curves intersect in Figures 17 and 20 (31), 
whereas the corresponding values for noncomplexed 
glutamyl species were obtained from the chemical 
shift vs. pD curves for glutathione at pD 7 where the 
amino group is protonated and the glutamyl carboxyl 
group is ionized. It can be seen from Table XVIII that 
the values are larger for zinc than for cadmium, in 
agreement with the relative order of stabilities of the 


corresponding glutamic acid complexes of these metal 
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ions (41). It is also evident that binding at the 

glutamyl end of glutathione affects the chemical shift 

of the Glu-COOH carbon more than that of the Glu-C,, 

carbon, which in turn experiences a larger effect 

than the Glu-C carbon. The order of these chemical 

shift effects is similar to that observed on protonation 

of the amino group in amino acids and peptides (27). 
The A values for the Gly-C, carbon listed in 

Table XVIII were obtained from the differences between 

the solid’ and dashed curves at pD 7 in Figures 17, 

20, and 22, and reflect the extent of binding of the 

glycyl carboxyl group of glutathione by these metal 

ions. It can be seen that the order of decreasing 

A values is pb*t > eat > ee the same as the relative 

order of the formation constants of the complexes 

of these metal ions with acetylglycine and the carboxyl 

group of polyglycine peptides (15,31). These carboxyl 

interactions are very weak as evidenced by the small 

A values, which are all less than 0.3 ppm. Binding of 

Zinc to the carboxyl groups of glutathione has not 

been previously detected in pH titration experiments 


(19=22)* presumably because the formation constants 


are so small. 
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Conclusions 
1 : 

C chemical shift data 
ions zinc, cadmium, lead and 
potential coordination sites 
high degree of specificity. 
bind to the sulfhydryl group 


cadmium bind to the glutamyl 


conditions used in this work. 
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indicate that the metal 
mercury bind to the 

of glutathione with a 
All four metal ions 
while only zinc and 
amino group under the 


These results demonstrate 


the potential of cmr as a technique for elucidating 


the binding of metal ions by 


should also be applicable to 


peptides. The technique 


Similarcostuditesefortother 


multidentate biological molecules. 
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CHAPTER VII 


EXPERIMENTAL 


A. CHEMICALS AND SOLUTIONS 


Reduced glutathione (Nutritional Biochemicals 
Corp. and Terochem Laboratories) was washed with a 
water-ethanol solution and dried at 110° before using. 
Reagent grade metal nitrate salts were used as 
received. 


The solutions were prepared in D.O or triply 


2 
distilled water under an atmosvhere of nitrogen to 

Minimize oxidation of the sulfhydryl group of gluta- 
thione. For the chemical shift vs. pD measurements, 


Ene pbrotL DAO solutions, 0.30 M in glutathione. and 


“ 
0.15 M in metal nitrate, was adjusted with a 40% 

KOD solution and concentrated nitric acid and samples 
of about 2 ml each were withdrawn at the appropriate 
pD values. 1,4-dioxane was added to the solutions 

as an internal reference at a concentration of about 
0.1 M. For making chemical shift vs. mole ratio 
measurements at constant pH, the requisite amount of 
metal nitrate was added to an HO solution of gluta- 
thione, the pH adjusted and a sample withdrawn. More 


metal nitrate was added to the solution and the same 


procedure followed for all samples at this pH. 


bp 
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B. POTENTIOMETRIC AND CMR MEASUREMENTS 


All pH and pD measurements were carried out as 
described in Chapter IV. 

The ide spectra were obtained using a Bruker 
HFX-90 spectrometer operating at a frequency of 22.63 
MHZ and equipped with a Nicolet 1085 computer. The 
Fourier transform mode was used with proton decoupling. 
When D,0 was the solvent, the deuterium resonance 
from the D.O was used for the heteronuclear lock signal. 
When HO was the solvent, the a resonance from hexa- 
fluorobenzene in a coaxial capillary was used for the 
lock. For each free induction decay signal, 8K 
data points were collected in the computer and 4K 
etna tone were carried out to achieve an adequate 
signal-to-noise ratio. The frequency range of the 
transformed spectra was 5000 Hz. Chemical shifts 
are reported in ppm relative to the resonance of 
internal 1,4-dioxane which is 67.4 ppm downfield from 
the TMS resonance. Positive chemical shifts correspond 
to greater shielding than in 1,4-dioxane. The chemical 
shift measurements are considered accurate to 
within 0.1 ppm. For all measurements the sample 


temperature was 32 = 2°. 
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